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f->^ ■ Abstract 

QQ ' We establish a long time soliton asymptotics for a nonlinear system of Maxwell equation coupled to 

f^ , a charged particle. The coupled system has a six dimensional manifold of soliton solutions. We show 

that in the large time approximation, any solution, with an initial state close to the solitary manifold, 

is a sum of a soliton and a dispersive wave which is a solution of the free Maxwell equation. It is 

^^ ' assumed that the charge density satisfies the Wiener condition. The proof is based on a development 

^ . of the general strategy introduced in the papers of Soffer and Weinstein, Buslaev and Perelman, and 

others: symplectic projection in Hilbert space onto the solitary manifold, modulation equations for 

the parameters of the projection, and decay of the transversal component. 
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1 Introduction 

Our contribution investigates the long time behavior of a single relativistic charge coupled to the Maxwell 
field. It is convenient to write the equations of motion in Hamiltonian form. The dynamical variables are 
then the position, q, of the particle, together with its momentum P, and the transverse vector potential, 
A, together with the canonically conjugate field Eg, which physically corresponds to transverse electric 
field. We refer to [371 Chapter 13] for details. In these variables the Hamiltonian function reads 



H{Es, A,q,P) = ]^j{\E,\'^ + I VAp) d^x+ l+[P- j p{x - q)A{x)d^x 



1/2 

(1.1) 



The phase space will be defined below through H < oo. The canonical equations of motion follow as 

(1.2) 



(1.3) 



s{x,t) = 


-AA(x,t) - 


-Il,{p{x-q{t))v{t)), A{x, 


t) = 


-Esix,t) 


rr(f) — 


P{t)- 


- p{x - q(t))A{x,t)(fx 


=■■ v{t), 


^v-) 


^+{P{t)- 


- / p{x - q{t))A{x,t)(fx] 
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3 
P{t) = \\, 


Jkit) [ p(x-q(t))VAk(xA) 


d^x. 





(1.4) 



k=l 

together with the transversality conditions 



V ■E,{x,t) = 0, V-A{x,t)=0. (1.5) 

Here and below all derivatives are understood in the sense of distributions. Es{x,t) is the projection of 
electric field onto the space of solenoidal (divergence-free) vector fields. In Fourier space the corresponding 
projector reads 

A a ■ k 

Usa = a t:2~^- 

p is the charge distribution of the particle, on which we will comment below. We use units such that the 
velocity of light c = 1, Eq = 1, and the mechanical mass of the charge m = 1. 
We note that, because of translational invariance, the total momentum 

ViEs, A,q,P)=P+ f Es{x) A (V A A{x))d^x (1.6) 



is conserved along sufficiently smooth trajectories of ()1.2p - (|1.5p . 
Let us write the system (|1.5p - (jl.4p as 

Y{t) = F{Y{t)), tGlR, (1.7) 

where Y{t) := {Es{x,t), A{x,t), q{t), P{t)). Below we always deal with column vectors but often write 
them as row vectors. The system (ll.2p - (ll.5p admits special solutions where the charge travels with 
constant velocity. In analogy with travelling solutions of nonlinear wave equation we call them solitons. 
Explicitly they are given by 



Ya,vit) = {Es,v{x - vt - a),A^{x - vt - a),vt + a,Py), P^ = v/^/l - v"^ + {p,A^) (1.8) 

for all a, w G IR with 1^1 < 1. The states Sa^v '■= ^a,t;(0) form the solitary manifold 

5:={Sa,^, :a,^;GlR^|^;| < 1}. (1.9) 



For general initial data one expects that for large times the solution splits up into two parts: one piece 
consists of a soliton with a definite velocity and the second piece are scattered fields escaping to infinity. 
In fact, this will be our main result. If the initial data are close to the solitary manifold, then we will 
prove that for large t 

{Es{x,t),A{x,t)) ~ {Es,v±{x - v±t - a±), Ay^{x - v±t - a±)) + W°{t)^±, t ^ ±00. (1.10) 

Here W^{t) is the dynamical group of the free wave equation (Equations (jl.2p with p = and (jl.Sp ). 
^± are the corresponding asymptotic scattered fields, and the remainder converges to zero in the global 
energy norm, i.e. in the norm of the space J- := //^(IR ) © Hl(Si ), see Section 2. For the particle 
trajectory we prove that 

q{t) ^^ v±, q{t)^v±t + a±, t ^ ±00. (l-H) 

The results are established under the following conditions on the charge distribution: p is a real valued 
function of the Sobolev class H^(JR ), compactly supported, and spherically symmetric, i.e. 

p,Vp,VVp£ L^iM-^), p{x) = iov \x\ > Rp, p{x) = pi{\x\). (1.12) 

An essential point of our asymptotic analysis is the Wiener condition 

p{k) = (2^)-3/2 f e'^'-'p{x)dx / for all A: G IR^ \ {0} . (1.13) 

The Wiener condition was noted already in the previous works [16\ [20l [2T\ [2^ [25} [26] . It expresses that 
all modes of the Maxwell field are coupled to the particle. 

There is no restriction on / \p{x)\dx. However if J p{x)dx 7^ 0, then the fields have a slow decay as 
|x|~2 at infinity. With our methods such a decay seems to be difficult to control and we have to impose 
the charge neutrality condition 

p{x)dx = 0. (1.14) 



(J1.14p is a technical condition. Physically one expects (jl.lOp to hold even without imposing charge 
neutrality and it is of interest to extend our proof in this direction. 

Let us briefly comment on earlier works. The first mathematical investigation is the contribution 
of Bambusi and Galgani [4j. They consider a non-relativistic kinetic energy for the charge and prove 
orbital stability of the solitons. The issue was taken up again in |16[ [271 I17j . where the kinetic energy 
is taken to be relativistic. In [16] the soliton-type asymptotics for the fields is established, under the 
Wiener condition, in local energy semi-norms centered at the particle's position q{t). The relaxation of 
acceleration, q{t) — > as i — > 00, is proved in [27] under the Wiener condition. Scattering behavior of 
the type (jl.lOp is established in [T7] for weak coupling, i.e. \\p\\l^ ^ 1. The case of a spinning charge is 
studied in [19] . 

A long time asymptotics as (ll.lOp appears also in nonlinear wave equations, like KDV and the U{1)- 
invariant nonlinear Schrodinger equation. Of course, in these equations there is no particle degrees 
of freedom, and S corresponds to the solitary wave solutions travelling at constant velocity. In fact 
our approach relies on and further develops general strategy introduced in [6l [71 [291 EO] • The strategy 
originates from the techniques of the paper [31] and their developments in [311 [32] [33] in the context 
of t/(l)-invariant Schrodinger equation. The approach uses the symplectic geometry methods for the 
Hamilton systems in Hilbert spaces and spectral theory of nonself adjoint operators. 

The invariant manifolds arise automatically for equations with a symmetry Lie group [5l [T3l [T3] . In 
particular, our system ()1.2l) - (ll.5p is invariant under translations in IR . The asymptotic stability of the 
solitary manifold is studied by a linearization of the dynamics ()1.7p . The linearization will be made on a 
special curve on the solitary manifold, S{t), which is the symplectic orthogonal projection of the solution. 
Then the linearized equation reads 

X{t) = A{t)X{t), t£M, (1.15) 



where the operator A{t) corresponds to the hnearization at the soliton S{t). Furthermore, consider the 
"frozen" hnearized equation (jl.lSp with A{ti) instead of A{t). The operator A{ti) has zero eigenvalue, 
and the frozen hnearized equation admits hnear in t secular solutions (see (j6.24p ). The existence of 
these runaway solutions prohibits the direct application of the Liapunov strategy and is responsible for 
the instability of the nonlinear dynamics along the manifold S. One of crucial observation is that the 
linearized equation is stable in the symplectic orthogonal complement to the tangent space Ts- The 
complement is invariant under the linearized dynamics since the linearized dynamics is Hamiltonian and 
leaves the symplectic structure invariant. 

Our proofs are based on a suitable development of the methods [61171129^130). Let us comment on the 
main steps. 

I. First, we construct the symplectic orthogonal projection S{t) = TlY{t) of the trajectory Y{t) onto 
the solitary manifold S. This means that S{t) E S, and the complement vector Z{t) := Y(t) — S{t) is 
symplectic orthogonal to the tangent space Ts(t) for every t G IR: 

Z{t)\Ts(^t), teJR. (1.16) 

So, we get the splitting Y(t) = S{t) + Z{t) and we linearize the dynamics in the transversal component 
Z{t) along the trajectory. 

The soliton component S{t) = Si,u\^u\ satisfies a modulation equation. Namely, in the parametrization 

^(t) = {c{t),v{t)) with c{t) := b{t) — / v{s)ds, we have 



Jo 

m = NMt),z{t)), \N,{m,zm<c\\z{t)f_f,, (i.i?) 

where || • ||_/3 stands for an appropriate weighted Sobolev norm. On the other hand, the transversal 
component satisfies the transversal equation 

Z{t) = A{t)Z{t) + N2{S{t), Z(t)), (1.18) 

where A{t) = ^5(4), and N2{S{t),Z(t)) is a nonlinear part: 

\\N2{s{t),zm\/,<c\\z{t)ef„ (1.19) 

where || • ||/3 is defined similarly to || • W-jS- Let us note that the bound ()1.19p is not a direct consequence 
of the linearization since the function S{t) generally is not a solution of (jl.7p . The modulation equation 
and the bound ()1.17p play a crucial role in the proof of (I1.19p . 

IL The linearized dynamics (jl.lSp is nonautonomous. First, let us fix t = ti in A(t) and consider the 
corresponding "frozen" linear autonomous equation with A{ti) instead of A(t). We prove the decay 

ll^(t)ll-/3<^|f|g^, teu (1.20) 

of the solutions X{t) to the frozen equation for any X{0) £ Zs^ where Si := S{ti), and Zs^ is the space 
of vectors X which are symplectic orthogonal to the tangent space Ts^ . Let us stress that the decay holds 
only for the solutions symplectic orthogonal to the tangent space. Basically, the reason of the decay is 
the fact that the spectrum of the generator A{ti) restricted to the space Zs^ is purely continuous. 

in. We combine the decay (|1.20p with the bound (|1.17p through the nonlinear equation (jl.lSp . This 
gives the time decay of the transversal component 

||Z,)||_, < £|», ..R, (Ol) 



if the norm ||Z(0)||/3 is sufficiently small. One of the main difficulties in proving the decay ()1.2ip is the 
non-autonomous character of the linear part of (jl.lSp . We deduce the decay from the equation (jl.lSp 
written in the "frozen" form 

Z{t) = A{ti)Z{t) + [A{t)-A{ti)]Z{t) + Ni{Sit),Z{t)), 0<t<ti, (1.22) 

with an arbitrary large ti > 0. 

IV. The decay ()1.2ip implies the soliton asymptotics (jl.lOp and (jl.lip by the known techniques of 
the scattering theory. 

Remarks 1.1 i) The asymptotic stability of the solitary manifold S is caused by the radiation of energy 
to infinity which appears as the local energy decay (11.2ip . 

ii) The asymptotics (jl.lOp can be interpreted as the collision of the incident soliton, with a trajectory 
V-t + a_, with an incident wave Wo{t)^^, which results in an outgoing soliton with a new trajectory 
^+^ + 0+, and a new outgoing wave Wo{t)'^+. It suggests to introduce the (nonlinear) scattering operator 

S: (7;_,a_,*_)^(i;+,a+,*+). (1.23) 

However, a correct definition of the operator is an open problem as well as the question on its asymptotic 
completeness (i.e. on its range). 

Remarks 1.2 i) The strategy [H [291 130] was developed also in [H [9l lOl El [Ml El ES |36] . Let us 

stress that all the papers contain several assumptions on the discrete and continuous spectrum of the 
linearized problem. In our case a complete investigation of the spectrum of the linearized problem is 
given under the Wiener condition, and there is no need in a priori spectral assumptions, 
ii) Note that the Wiener condition is indispensable for our proof of the decay (ll.20p . and only in one 
point: in the proof of Lemma 115.21 In all other places we use only the fact that the coupling function 
p{x) is not identically zero. The other assumptions on p can be weakened: the spherical symmetry is 
not necessary, and one can assume also that p belongs to a weighted Sobolev space rather than has a 
compact support. 

Our paper is organized as follows. In Section 2, we formulate the main result. In Section 3, we introduce 
the symplectic projection onto the solitary manifold. The linearized equation is defined in section 4 and 
studied in the next one. In Section 6, we split the dynamics in two components: along the solitary 
manifold, and in transversal directions, and we justify the estimate (|1.17p concerning the tangential 
component. The time decay of the transversal component is established in sections 7-10 under an 
assumption on the time decay of the linearized dynamics. In Section 11, we prove the main result. 
Sections 12 - 18 fill out the gap concerning the time decay of the linearized dynamics. In Appendices A 
to C we collect some routine calculations. 

Acknowledgments The authors thank V.Buslaev for numerous lectures on his results and fruitful 
discussions. 

2 Main Results 

2.1 Existence of Dynamics 

Let us introduce a phase space for the system (II. Sp to (II. 4p and state the existence of dynamics. Set H^ = 
L2(IR3,IR3), ifi is the closure of C^(IR^IR^) with respect to the norm P||i = |V^| = ||V^||^2(]e^3 j^3). 

Let Hg, Hi be the subspaces constituted by solenoidal vector fields, namely the closure in H^ , H^ 
respectively of C^ vector fields with vanishing divergence. Define the phase space 

£ = H^,® H] © IR3 © IR^ 



where the norm of y = {Eg, A, q, P) is 

||y||£ = 1^,1 + Pill + |g| + |P|. 

Let us define the corresponding space for fields alone: 

T = H'',(B Hi \\{Es, A)y = \Es\ + P||i. 

We write the Cauchy problem for the system (jl.5p - (|1.4p as 

Y = F{Y{t)), tGM; Y{0) = Y^. (2.1) 

Proposition 2.1 Let (C) holds, let Y° = {E^ , A^ , q^ , P°) G £. Then 

i) There exists a unique solution Y{t) £ C{1R,£) to the Cauchy problem 112. 1\) . 

ii) The energy and the total momentum are conserved, 

H{Y{t)) = H{Y°), V{Y{t))=V{Y^), t G M. 

Hi) The estimate holds, 

\q{t)\ <v <1, t£TR. (2.2) 

For proof see |16| . 



2.2 Soliton Solutions and the Main Result 

We will investigate the long-time behavior of finite energy solutions to (|1.5p - (|1.4p . The set of asymptotic 
solutions corresponds to the charge travelling with a uniform velocity, v. Up to translation they are of 
the form 

Es{x,t) = Eg^vix — vt), A{x,t) = Ay{x — vt), q{t) = vt, Py = p^, + / pAyd x, (2-3) 

where v £ IR^, \v\ < 1; Es,v = HsE'", A^ = UsA" , and E"" , A", p^ are given by 

E^(x) = -VMx) +V-V A^ix), A^ix) = v^y{x); 

S{x)-^[ P(^)^'^ „ 

YV\X) — / I / Nil' \ I ' "'^ 



4vr y \l{y-x)\\ + {y-x)±_y ^/l 



v^ 



Here 7 = l/Vl — v'^ and x = x\\+ xj_, where xii \\v and x_\_^-V for x G IR . Below we call the solutions 
of type (|2.3p "solitons" . For simplisity of notations let us omit the index s and write E^ instead of Eg^v 
The stationary equations for Ey, Ay,Py read 

Ey{x)=vVAy{x), (2.4) 

vVEy{x) = ^Ay{x)+Iis{p{x)v), (2.5) 

Pv - I p{x)Ay{x)d^x 

(2.6) 



1/2 ' 



l + (n- f pix)Ay{x)d^x^^ 
0= p{x)vVAy{x)d^x. (2.7) 



To state our main result we have to introduce the following weighted Sobolev spaces. Let H^ ^, H] ^ be 
the subspaces of -ff^ , respectively Hi consisting of all the fields E, resp. A with the finite norms 

\\E\\o,a = 1(1 + \x\rEl P||i,„ = 1(1 + \x\rVA\. 

Since H] C L6(]R3^]R3), we have H] C F.^, a < -1. Set 

£a = Hl^®Hl^®M^®M^, \\Y\\o, = \\Esh,a + \\A\\i,a + \q\ + \P\, Y&E-oc- 
For the fields we set 

^cx = Hg^ © Hg^^, \\{Es,A)\\a = ||£'s||o,o + ||^||l,o- 

By (|2.3p and the neutrality condition (J1.14p we obtain that E^ G H^,a for a < 3/2, A^ G H]^^ for a < 1/2 
and thus, 

Sa,v G £a-, for a < 1/2, where Sa,v '■= {E^{x — a)^A^{x — a),a,Pv). (2-8) 

The main result of our paper is the following theorem. 

Theorem 2.2 Let the condition lil.l^} . Wiener condition il.lcl]} . and the neutrality condition ^-14^ 
hold, let Y^ £ £p with P = l + 6, 0<6< 1/2. Suppose that Y^ is sufficiently close to the solitary 
manifold: 

Y'' = Sa,,vo + ZQ, dp:=\\Z4p^l. (2.9) 

LetY{t) G C{1R,£) be the solution to the Cauchy problem h2.1\) . Then the asymptotics hold fort -^ ±00, 

q{t) = ^± + Oi\t\-^~^), q{t)=v±t + a± + 0{\t\-^^y, (2.10) 

{Es{x, t),A{x, t)) = {E^^ {x - v±t - a±), A^^ {x - v±t - a±)) + W°{t)^± + r±{x, t) (2.11) 

with 

\\r±{t)\y = 0{\t\-'). (2.12) 

It suffices to prove the asymptotics (12. IIP . (12.101) for t -^ +cxd since the system (ll.5l) - (ll.4p is time 
reversible. 

3 Symplectic Projection 

3.1 Symplectic Structure 

From now and later on we omit the subindex s in hamiltonian variables and write Y = (E, A, q, P). The 
system (jl.Sp to (jl.4p reads as the Hamiltonian system 

/ 1 \ 
Y = JVn{Y), J:= ~o^ 1 ,Y = {E,A,q,P)e£, (3.1) 

V 0-10/ 

where DTC is the Frechet derivative of the Hamilton functional (jl.ip . Let us identify the tangent space 
to £^, at every point, with £. Consider the symplectic form il defined on £ by 



n= dE{x) A dA{x) dx + dqA dP, 



I.e. 



n{YuY2) = j{Ai •E2-A2- Ei)dx + qi-P2-q2-Pi (3.2) 

for Yfc = (E'fc, A}^,qk,Pk) G iS, A; = 1, 2 if the integral converges. 

Definition 3.1 i) Yi \ I2 means that Yi ^ £ is symplectic orthogonal to Y2 G £, i.e. ^{Yi,Y2) = 0. 

a) A projection operator P : £ ^ £ is called symplectic orthogonal if Yi f I2 for Yi G KerP and 
Y2 G ImP. 



3.2 Symplectic Projection onto Solitary Manifold 

Definition 3.2 A soliton state is S{a) := {Ey{x — b),Ay{x — b),b,Pv), where a := (6,^1) with 6, f G IR 
and \v\ < 1. 

Obviously, the soliton solution admits the representation S{a{t)), where 

a{t) = {b{t),v{t)) = {vt + a, v). (3.3) 



Definition 3.3 A solitary manifold is the set S := {S{b,v) : b G IR^, \v\ < 1}. 

Let us consider the tangent space 7^(o-)5 to the manifold 5 at a point S{a). The vectors tj := da.S{a)^ 
where d^ := dh, and So- _,_3 := d^. with j = 1, 2, 3, form a basis in 7^5. In detail, 



Tj = Tj{v) := db^S{a) = ( -djE^{y) , -djA^{y) , Cj , ) 



j = 1,2,3, (3.4) 



where y := x — b is the "moving frame coordinate", ei = (1,0,0) etc. Let us stress that the functions Tj 
will be considered always as the functions of y, not of x. 

Note that the symplectic form Q{ti{v) , Tj {v)) is well defined by the neutrality condition (J1.14p . 

Lemma 3.4 The matrix with the elements 0,{ti{v),tj{v)) is non- degenerate for \v\ < 1. 

The proof is made by a straightforward computation, see Appendix A. 

Now we show that in a small neighborhood of the soliton manifold S a "symplectic orthogonal 
projection" onto S is well-defined. Let us introduce the translations 

Ta : {ip{-),7ri-),q,p) ^ {ip{- -a),vr(- -a),q + a,p), a G IR^. 

Note that the manifold S is invariant with respect to the translations. 

Definition 3.5 Let us denote by v{Y) := P/Vl + P^ where P G IR is the last component of the vector 
Y. 

Lemma 3.6 Let [rW\) hold, -3/2 < a < 1/2 andv <1. Then 

i) there exists a neighborhood Oa{S) of S in Ea and a map H : Oa{S) -^ S such that H is uniformly 

continuous on Oa{S) n {Y G £a ■ v{Y) < v} in the metric of £a, 

UY = Y for Y eS, and Y - S\ TsS, where S = UY. (3.5) 

ii) Oa{S) is invariant with respect to the translations Ta, and 

UTaY = TJ1Y, for Y e Oo,{S) and a G IR^. (3.6) 

Hi) For any v < 1 there exists a v <1 s.t. |f (ny)| < v when \v{Y)\ < v. 

iv) For any v < 1 there exists an ra{v) > s.t. S{a) + Z £ Oa{S) if \v{S{a))\ < v and \\Z\\a < ra{v). 

The proof is similar to that of Lemma 3.4 in |20j . 

We will call II the symplectic orthogonal projection onto S. 

Corollary 3.7 The condition ^2. 9\) implies that Yq = S + Zq where S = S{crQ) = HYq, and 

\\Zoh<^l. (3.7) 



4 Linearization on the Solitary Manifold 

Let us consider a solution to the system ()1.5p - (|1.4p . and split it as the sum 

Y{t) = S{a{t)) + Z{t), (4.1) 

where a{t) = {b{t),v{t)) € IR x {\v\ < 1} is an arbitrary smooth function of t G IR. In detail, denote 
Y = {E,A,q,P) and Z = (e, a, r, vr). Then ()4.ip means that 

r,{+\ — hd^ J- 7-(+\ 

(4.2) 



E{x,t) = E^^t){x-b{t))+eix-bit),t), q{t) = b{t)+r{t) 
A{x,t) = A^^t){x-b{t))+aix-b{t),t), Pit) = P^(i)+7r(t) 



Let us substitute (|4.2p to (|1.5p - (|1.4p and linearize the equations in Z. Later we will choose S{a{t)) = 
TIY{t), i.e. Z{t) is symplectic orthogonal to Tsr^u\\S. However, this orthogonality condition is not 
needed for the formal process of linearization. The orthogonality condition will be important in Section 
6, where we derive "modulation equations" for the parameters cr{t). 

Let us proceed to linearization. Setting y = x — b{t) which is the "moving frame coordinate", we 
obtain from ([421) and ((OJ-dLil) that 

E = v-V^E^^t){y)-b-VE^^t){y) + e{y,t)-b-Ve{y,t) 

= -A{A,{y)+a{y,t))-UMy-r)q), (4.3) 

A = v- V^^^(i) (y) - 6 • VA^^t) (y) + d(y, t)-b- Va{y, t) = -E^^t) iv) - e(y, t), (4.4) 

. ^ ; . ^ Pyjt) + TT - {pjy - r), A^(t) (y) + a{y, t)) 

^" ^'■"(l + (P.(t)+^-(p(y-r),^(,)(y)+G(y,t)))2)i/2' ^ ' > 

P = v- V„P„(t) +^ = Y^ Qkipiy - r), V(A^(i)(y) + a{y, t))k). (4.6) 

k 

Step i) First we linearize the equation (|4.5p . Note that p{y — r) = p{y) — r ■ Vp{y) + N2{r), where 

\N2{r)\f3 < Cp{ry, (4.7) 

uniformly in \r\ <r for any fixed r. Here (3 is the constant defined in Theorem 12.21 Then (let us write v 
instead of v{t) and omit the other arguments for simplicity) 

(p(y - r), A^ + a) = {p, A^) + {p, a) - {r ■ Vp, A^) + N^ = {p, A„) + (p, a) + N^, (4.8) 

where N2{r, a) = — (r • Vp, a) + {N2,Ay + a). Here we use the equality (r • Vp, Ay) = which holds, since 
A^ is even and Vp is odd. Further, since P^ — (p, A^) = Pv by (j2.3p . we get Pj, + vr — {p{y — r),Ay + a) = 
Py + TT - (p. Ay) - {p,a) - N2 = py + TT - (p, a) - N2 = Pv + s, where s := vr - (p, a) - N2. Applying 
Taylor expansion we obtain 

(1 + {p. + s)')~'/' = \ - ^- ,^3 +n:,= .^^ \. ^ ,, - -^ + N',, 
{i+piyi^ [i+pif^ [i+piy/^ i + Pt, 

since Pv/i'i- +p1Y = v. Finally, 

Py + Ti- {p{y - r),Ay + a) Pv + s 



(1 + (P, + vr - {p{y - r),Av{y) + a))2)i/2 (1 + {py + s^/^ 



where u = u^ = {1 — v'^Y''^ = (1 +p^) ^'^. Insert the expression for s, then the equation (j4.5p becomes 

f = v-b + B^{7T-{p,a))+Ns, (4.9) 

where By := z^(£' — f (8) t>), and 

\N,{Z)\<C{v)\Z\lf, (4.10) 

uniformly in |t)| < u < 1. 

5iep ii) Next we hnearize the equation (14. 3p . By (14.81) and (|4.9p we obtain /9(y — r)q = pv + pBy{-ir — 
{p, a)) — r ■ V pv + N[. Substitute to the equation (j4.3p and take (|2.4p into account, then we get 



e = b-Ve- Aa + {b-v)- VE^ - v ■ VyEy - n^(pS„(7r - {p,a)) - r ■ V pv) + iVi, 

where for A^i the same bound (I4.10p holds. 

Step in) Further, by (j2.5p the equation ()4.4p becomes 

a = -e + b ■ V a + {b - v) ■ VAy - v ■ V^Ay. 

Step iv) Let us proceed to the equation (j4.6p . We have 

P = v. VyPy + n= {p{y -r),q- V{Ay + a)) 

= {p-r-Vp + N2,{v + By{7r - {p, a)) + iVg) • V(A„ + a)) 

= {p-r-Vp + N2,v- V{Ay + a) + By{TT - {p, a)) ■ VAy + N^) 

= {p,v- V{Ay + a)) + {p, By{TT - {p, a)) ■ VAy) -{r-Vp,v- VAy) + N4 

= {p,v ■ Va) - {r ■\7p,v- VAy) + N4, 

since (/>, f • VAy) = and {p,By{i: — {p,a)) • VAy) = 0. Finally, the equation becomes 

jr = {p^v Va) - V • VyBa - {r • Vp, v ■ VAy) + 7V4, 

where for Ni{v, Z) the estimate similar to (j4.10p holds. 
We can write the equations (IT9l) . (|4.11l) - ()4.13p as 

Z{t)=A{t)Z{t)+T{t) + N{t), tGlR. 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



Here the operator A{t) depends on a{t) = {b{t),v{t)). We will use the parameters v = v{t) and w := b{t). 
Then A{t) can be written in the form 



/-\ 



Ait) 



a 
r 



Ay, 



\^ J 



( e\ 

a 

r 

V vr/ 



f wV -A + Us{pBy{p,-)) Usi-Vpv) -UsipBy) \ f e\ 

-1 wV 

-By{p,-) By 

\ {p,v-V-) -{-Vp.vVAy) 

Furthermore, T[t) and N{t) in (j4.14p stand for 



a 
r 



J 



(4.15) 



v^y 



T{t) = Ty 



/ {W-V)- VEy - V ■ VyEy \ 
{W — V) ■ VAy — V ■ VyAy 

V — W 

\ -V ■ VyPy ) 



N{t) = N{a,Z) 



/ N,{v,Z) \ 

N2{Z) 
Ns{v,Z) 
\ N,iv,Z) J 



(4.16) 



where v = v{t), w = w{t), a = a{t) = {b{t),v{t)), and Z = Z{t). The estimates ()4.7p and (|4.10p . which 
holds as weh for A^i and A'^4 imply that 

\\N{a,Z)y<C{v,Q)\\Z\tp (4.17) 

uniformly in a and ||Z||_^ < r„^({i) for any fixed v < 1 and \v\ < v. 

Remarks 4.1 i) The term A{t)Z{t) in the right hand side of the equation (j4.14p is linear in Z(t), and 
N{t) is a high order term in Z{t). The term T{t) vanishes if S{a{t)) is a soliton solution since in this 
case V = and w = b = v. Otherwise T{t) is a zero order term which does not vanish though S{(T{t)) 
belongs to the solitary manifold, 
ii) Formulas ()3.4p and (j4.16p imply: 

3 

T{t) = -Y^[{w-v)iTi + viTi+3] (4.18) 

1=1 

and hence T{t) £ T^^^^^^^iS, i E IR. 

5 The Linearized Equation 

Here we study some properties of the operator (j4.15p . First, let us compute the action of Ay^^u on the 
tangent vectors Tj to the solitary manifold S. 

Lemma 5.1 The operator A^^^ acts on the tangent vectors Tj{v) to the solitary manifold as follows, 

Av,w[Tjiv)] = {w -V) ■ VTj{v), Ay^w[Tj+3{v)] = {w - v) ■ Vtj+3{v) + Tj{v), j = 1, 2, 3. (5.1) 

Proof We differentiate the stationary equations (|2.4p to (j2.7p and obtain 

djEy = V ■ VdjA^, V ■ VdjEy = AdjA^ + Us{djpv), (5.2) 

dy^Ey = djAy + V ■ Vdy^Ay, OjEy + V ' VOy^Ey = AOy^Ay + UsipCj), (5.3) 

Cj = Bydy^Py = ByOy^ {Py " (/), Ay)), = (p, Bj " VAy) + {p,V ' VOy^Ay). (5.4) 

Let us apply Ay^^ to tj, j = 1,2,3, i.e. to {—djEy,—djAy,ej,0). The first component equals 

w ■ V{-djEy) - A{-djAy) + Us{pBy{p, " ^j' ^ ) ) + Us{djpv) = {w-v)- V{-djEy) (5.5) 

by (|5.2p and since {p,djAy) = 0. The second component equals 

djEy - w ■ V{djAy) = {w-v)- V{-djAy) (5.6) 

by (|5.2p . The third component equals 

-By{p,-d,Ay)=0, (5.7) 

and the last is 

{P,V V{-djAy)) - {p,By{p,-djAy) ' V Ay) " {djP,V V Ay) = 0. (5.8) 

The formulas (j5.5p to (j5.8p mean j4^,^^„[rj] = {w — v) ■ Vtj, j = 1, 2, 3. 

Now let us apply Ay^^ to tj_|_3, j = 1, 2, 3, i.e. to {dy.Ey,dy^Ay, 0, dy.Py). The first component equals 

W ■ Vdy^Ey - A{dy^Ay) + H^ipBy {p , Oy^ Ay)) 
= {W -V) ■ Vdy^Ey - djEy + UsipCj) + Us{pBy{p, dy^Ay)) " Us{pBydy^Py) 
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= {w-v)-Vd^^E.,-djE^ (5.9) 

by (|5.3p and (j5.4p . The second component equals 

-dy.E^ + w Vdy. A^ = {w-v)- V5„^. ^„ - djA^ (5.10) 

by (|5.3p . The third component equals 

- Bv{p, dy.A^) + BydvjPv = Bydv^pjj = ej (5-11) 

by (|5.4p . and the last is 

{p, V • V9^^. A„) - {p, B^{p, d^^A^) • VA„) + {p, (Bydy^Py) • VA) 

= {P,V- Vdy^Ay) - {p, By{dy^Py - By^Pv) ■ V Ay) + (/), {ByBy^Py) ■ V Ay) 

= {P,V- Vdy^Ay) + {p, djAy) = (5.12) 

by (IS3D- The formulas ([53]) to ([^1^ mean Ay^yj[Tj+3] = {w - v) ■ Vtj+s + r^, j = 1, 2, 3. D 

Consider the linear equation 

X{t)=Ay,y,X{t), telR (5.13) 

with an arbitrary fixed v such that |w| < 1 and w G IR . 

Lemma 5.2 i) For any v, \v\ < 1, w ^ IR the equation 115. 13\) formally can be written as the Hamiltonian 
system (of. i TO]] ). 

X{t) = JDHy^UXit)), teJR, (5.14) 

where DTCy is the Frechet derivative of the Hamilton functional 

'Hy^w{X) = - |e|2 + |Vap dy + a{w ■ V)edy + -{By{p, a)) ■ {p, a) 

+ -vr • ByTT + (r • Vpv, a) - {pByir, a) + -{r ■ Vp, v ■ (VAy)), X = (e, a, r, vr) G <S. (5.15) 

ii) Energy conservation law holds for the solutions X{t) G C^(iR,£), 

ny,yj{X{t)) = const, t G IR. (5.16) 

Hi) The skew-symmetry relation holds, 

n{Ay^y,x,,X2) = -n{x^,Ay^^X2), x^x^es. (5.17) 

We will apply Lemma 15.21 mainly to the operator Ay^y corresponding to w = v. In that case the 
linearized equation has the following additional essential features. 

Lemma 5.3 Let us assume that w = v and \v\ < 1. Then 

i) The tangent vectors Tj{v) with j = 1, 2, 3 are eigenvectors, and Tj^3{v) are root vectors of the operator 

Ay^y, corresponding to zero eigenvalue, i.e. 

AvA-^jiv)] = 0, AyATj+3{v)] = Tj{v), j = 1,2,3. (5.18) 

ii) The Hamilton function ( t5. J5|) is positive definite, 

Hy^X) > 0. (5.19) 

11 



Proof The first statement follows from (j5.ip . To prove the second statement let us rewrite 
n^A^) = 2 / |eP + |Vap dy+ a{v ■V)edy + -{B.,{p,a)) ■ {p,a) 

+-7r • B^TT + (r • Vpv, a) - {pB^n, a) + -{r ■ Vp, v ■ (VA^)), X = (e, a, r, vr) G £. 
= liBvi^r - {p, a))) ■ (tt - {p, a)) + ^{{e, e) + {{v ■ V)a, {v ■ V)a) - (e, {v ■ V)a)) 

+^(((-A + {v ■ Vf)a, a) + {{r ■ V)pv, a) + ((r • V)p, v ■ (V^)))- 

Here the first line is clearly nonnegative, since B^ is nonnegative definite. The last line in Fourier space 
by (|A.4p equals 

1 / ((A:^ - ikvfM' - mr)p{v • I) + ^^^) dk. 
The integrand is nonnegative, since 

'iramrmv-Tm e IR and |Re|i(tr)«„ .I)l| < |(tr)||/,|M|a|. 

This completes the proof. D 

Remark 5.4 Lemma 15.31 ii) together with energy conservation (j5.16p will imply the analyticity of the 
resolvent {A^^y — \)~^ for ReA > 0, see Section 14. 

Remark 5.5 For a soliton solution of the system (ll.5p - (ll.4p we have h = v, v = Q^ and hence T{t) = 0. 
Thus, the equation (j5.13p is the linearization of the system (jl.5P " (|1.4p on a soliton solution. In fact, we 
do not linearize ()1.5p ~ ()1.4p on a soliton solution, but on a trajectory 5(cj(t)) with a{t) being nonlinear in 
t. We will show later that T{t) is quadratic in Z{t) if we choose S{a{t)) to be the symplectic orthogonal 
projection of Y{t). Then (I5.13P is again the linearization of (|1.5p - (|1.4p . 

6 Symplectic Decomposition of the Dynamics 

Here we decompose the dynamics in two components: along the manifold S and in transversal directions. 
The equation ()4.14p is obtained without any assumption on a{t) in ()4.ip . We are going to choose 
S{a{t)) := ny(t), but then we need to know that 

Y{t) e 0„(5), t e IR, (6.1) 

with some Oa{S) defined in Lemma [3.5[ It is true for t = by our main assumption (12. 9p with sufficiently 
small d,3 > 0. Then S{a{0)) = nY{0) and Z(0) = Y{0) - S{a{0)) are weU defined. We will prove below 
that ()6.ip holds with a = —(3 if d/3 is sufficiently small. First, the a priori estimate (12. 2p together with 
Lemma 13.61 iii) imply that ny(t) = S{a{t)) with a{t) = (b{t),v{t)), and 

\v{t)\<v<l, telR (6.2) 

if Y{t) £ 0^p{S). Denote by r_p{v) the positive number from Lemma 13.61 iv) which corresponds to 
a = -13. Then S{a) + Z £ 0-f}{S) if cj = {b,v) with \v\ < v and \\Z\\_,3 < r_p{v). Note that ([22]) implies 
||Z(0)||_^ < r_p{v) if dp is sufficiently smah. Therefore, S{a{t)) = HY{t) and Z{t) = Y{t) - S{a{t)) 
are well defined for t > so small that ||Z(t)||_^ < r-j^{v). This is formalized by the following standard 
definition. 
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Definition 6.1 t^ is the "exit time", 

t^ = sup{t>0:\\Z{s)\\-f3<r^f3{v), < s <t}, Z{s) = Y{s) - S{a{s)). (6.3) 

One of our main goals is to prove tiiat t^, = oo if df^ is sufficiently small. This would follow if we show 
that 

||Z(t)||_^<r_^(i))/2, 0<t<t,. (6.4) 

Note that 

\r{t)\<r:=r_fs{v), < t < t^. (6.5) 

Now N{t) in (|4.14p satisfies, by (|4.17p . the following estimate, 

\\N{t)y<Cp{v)\\Zm-(3, 0<t<t,. (6.6) 

6.1 Longitudinal Dynamics: Modulation Equations 

From now on we fix the decomposition y{t) = S{a{t)) + Z{t) for < t < t,, by setting S{a{t)) = ny(i) 
which is equivalent to the symplectic orthogonality condition of type (j3.5p , 



Z(t) t Ts(,(i))cS, 0<t<t,. (6.7) 

This allows us to simplify drastically the asymptotic analysis of the dynamical equations (j4.14p for the 
transversal component Z{t). As the first step, we derive the longitudinal dynamics, i.e. the "modulation 
equations" for the parameters (y{t). Let us derive a system of ordinary differential equations for the 
vector (T{t). For this purpose, let us write (|6.7p in the form 

J](Z(t),r,(t))=0, j = l,...,6, 0<t<t*, (6.8) 

where the vectors Tj{t) = Tj{a{t)) span the tangent space Tg(^„(^f^)S. Note that a{t) = {b{t),v{t)), where 

\v{t)\<v<l, 0<t<t,, (6.9) 

by Lemma [3^ iii). It would be convenient for us to use some other parameters (c, v) instead of o" = (6, v), 

where c{t) = b{t) — / v{T)dT and 
Jo 

c{t) = b{t)-v{t) = w{t)-v{t), 0<t<t*. (6.10) 

We do not need an explicit form of the equations for (c, v) but the following statement. 

Lemma 6.2 (ci. [20j, Lemma 6.2) Let Y{t) be a solution to the Cauchy problem V2.1\) . and 1^. j| ), \6.8\) 
hold. Then {c{t),v{t)) satisfies the equation 

(S) )=-^(^(*)'^W)' 0<t<t„ (6.11) 

where 

M{a,Z) = 0{\\Z\t^) (6.12) 

uniformly in a ^ 5]('u) := {{b,v) : \v\ < v}. 
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Proof We differentiate tlie orthogonality conditions (j6.8p in t, and obtain 

o = n{z,Tj) + n{z,fj) = n{Az + T + N,Tj) + n{z,fj), o<t<u. (6.13) 

First, let us compute the principal (i.e. non-vanishing at Z = 0) term Q{T,tj). For j = 1,2,3 one has 
by dHHD, ([Alii) 

^{T, Tj) = - ^(QO(Ti, Tj) + Vin{Tl+3, Tj)) = ^ n{Tj,Tl+3)vi = ^ nj^Vl, 
I I I 

where the matrix fi"*" is defined by (1A.13|1 . Similarly, 

^{T, Tj+s) = - ^(Qfi(r/, Tj+s) + Vin{Ti+s,rj+3)) = ^ ^(^j+3, Tl)ci = - ^ (7^Q. 

I I I 

As the result, we have by ()A.14p . 

in the vector form. 

Second, let us compute i}{AZ,Tj). The skew-symmetry (j5.17p implies that Q{AZ,tj) = —Q{Z,Atj). 
Then for j = 1, 2, 3, we have by (15. ip . 



n{AZ, Tj) = -n{Z, c ■ Vtj), (6.15) 

and similarly, 

n{Az, Tj+s) = -o(z, c • Vtj+3 + Tj) = -o(z, c • Vtj+s) - n{z, Tj) = -n{z, t ■ Vt^+s), (e.ie) 

since il.{Z,Tj) = 0. 

Finally, let us compute the last term Q{Z, fj). For j = 1, . . . , 6 one has fj = b-VbTj + v-VyTj = v-VyTj 
since the vectors tj do not depend on h according to ()3.4p . Hence, 

^{Z,fj) = 9.{Z,v-VyTj). (6.17) 

As the result, by ()6.14p - ()6.17p . the equation ()6.13p becomes 

= n{v) (''A+Mq{cj,Z)(''A+ AAo(a, Z), (6.18) 

where the matrix A1o(o",Z) = OdlZH./?), and A/'o(o-, Z) = 0(11^11^^) uniformly in cr € T.{v) and ||Z||_/3 < 
r_^(i}). Then, since rt{v) is invertible by Lemma 13.41 and ||Z||_^ is small, we can resolve ()6.18p with 
respect to the derivatives and obtain equations (j6.1ip with J\f = 0{\\Z\iin) uniformly in a E 5](i;). D 

Remark 6.3 The equations (16. lip . ()6.12p imply that the soliton parameters c(t) and v{t) are adiabatic 
invariants (see [2]). 

6.2 Decay for the Transversal Dynamics 

In Section 11 we will show that our main Theorem 12.21 can be derived from the following time decay of 
the transversal component Z{t): 

Proposition 6.4 Let all conditions of Theorem \2.2\ hold. Then t^ = oo, and 

C{p,v,dfj) 
-"- {l + \t\)P ' 
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(6.19) 



We will derive (j6.19p in Sections 8-11 from our equation (j4.14p for the transversal component Z(t). This 
equation can be specified using Lemma |6.2[ Indeed, the lemma implies that 

||r(t)||^<C7(^))||Z(t)||2_^, 0<t<t„ (6.20) 

by (|4.16p since w — v = c. Thus (j4.14p becomes the equation 

Z{t) = A{t)Z{t)+N{t), 0<t<U, (6.21) 

where A{t) = ^„(t),«,(t), and N{t) := T{t) + N{t) satisfies the estimate 

\\N{t)\\^<C\\Z{t)f_f„ 0<t<U. (6.22) 

In all remaining part of our paper we will analyze mainly the basic equation (I6.21|) to establish the 
decay ()6.19p . We are going to derive the decay using the bound ()6.22p and the orthogonality condition 
dSH). 

Let us comment on two main difficulties in proving (j6.19p . The difficulties are common for the 
problems studied in [Jl [10]. First, the linear part of the equation is non- autonomous, hence we cannot 
apply directly known methods of scattering theory. Similarly to the approach of [3 |T0] , we reduce the 
problem to the analysis of the frozen linear equation, 

X{t) = AiX{t), t G IR, (6.23) 

where Ai is the operator A^-^^^y^ defined in (14.151) with vi = v{ti) and a fixed ti £ [0, t^,). Then we estimate 
the error by the method of major ants. 

Second, even for the frozen equation (j6.23p . the decay of type (|6.19p for all solutions does not hold 
without the orthogonality condition of type (j6.7p . Namely, by (jS.lSp the equation (j6.23p admits the 
secular solutions 

3 3 

X{t) = J^C,r,(^;i) + J]D,[r,(t;i)t + r,+3(«i)] (6.24) 

1 1 

which arise also by differentiation of the soliton (jl.Sp in the parameters a and vi in the moving coordinate 
y = X — fit. Hence, we have to take into account the orthogonality condition ()6.7p in order to avoid 
the secular solutions. For this purpose we will apply the corresponding symplectic orthogonal projection 
which kills the "runaway solutions" (I6.24p . 

Remark 6.5 The solution (I6.24p lies in the tangent space Tgf„^\S with ai = (6i,fi) (for an arbitrary 
hi G IR) that suggests an unstable character of the nonlinear dynamics along the solitary manifold (cf. 
Remark 14.11 iii)). 

Definition 6.6 i) Denote by 11^,, \v\ < 1, the symplectic orthogonal projection of £ onto the tangent 
space Ts(^^^S, and Py = I — 11^,. 

ii) Denote by Z^ = Py£ the space symplectic orthogonal to Tstfj\S with a = (6, v) (for an arbitrary 

6 sir;. 

Note that by the linearity, 

n,Z = ^n,i(7;)Tj(f)f^(r;(f),Z), Ze£, (6.25) 

with some smooth coefficients Y\.ji{v). Hence, the projector 11^, in the variable y = x — b, does not depend 
on b, and this explains the choice of the subindex in n„ and P^ . 
Now we have the symplectic orthogonal decomposition 

£l3 = '^S{a)'S + 2^, a = {b,v), (6.26) 

and the symplectic orthogonality ()6.7p can be written in the following equivalent forms, 

n,(t)Z(t) = 0, Py^t)Z{t) = Z{t), 0<t<t,. (6.27) 
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Remark 6.7 The tangent space Tg/^\S is invariant under the operator j4^„ by Lemma [5.31 i). hence the 
space Z^ is also invariant by (j5.17p : A^^^Z S Z^ for sufficiently smooth Z £ Z^j. 

In the Sections 12-18 we prove the fohowing proposition which is one of the main ingredients for 
proving (16.19p . Let us consider the Cauchy problem for the equation (|6.23p with Ai = Ay^^^^-^ for a fixed 
vi, \vi\ < 1. Recall that the parameter f3 = l + 6,0<6< 1/2 is also fixed. 

Proposition 6.8 Let the Wiener condition /il.l3\) and the neutrality condition (TTT^ hold, \vi\ <v<l, 

and Xq £ £. Then 

i) The equation I16.23\} . with Ai = Ay-^^^-^, admits the unique solution c^^Xq := X{t) G C(\R,£) with the 

initial condition X(0) = Xq. 

a) For Xq G Zy-^ n £j3, the solution X{t) has the following decay, 

||e^^%||~/3 < ^^^^ll^oll^, teM. (6.28) 

Remark 6.9 The decay is provided by two fundamental facts which we will establish below: 

i) the null root space of the generator Ai coincides with the 7^(o-^)5 with cji = {bi,vi) (for an arbitrary 

6i G M), and 

ii) the spectrum of Ai in the space Z^-^ is absolutely continuous. 

7 Frozen Form of Transversal Dynamics 

Now let us fix an arbitrary ti G [0, t*), and rewrite the equation ()6.2ip in a "frozen form" 

Z{t) = AiZ{t) + iA{t)-Ai)Z{t) + Nit), 0<t<t^, (7.1) 

where Ai = ^^(tj)^„(j^) and A{t) - Ai = 

f [w-vi]-V UMBv-By,){p,-)) Us{-Vpiv-vi)) -Usip{By - By,)-) \ 

[w-vi]-V 

-{By-By,){p,-) By-By, 

\ {p,{v-vi)V-) -{■Vp,{vVAy-viVAy,)) / 

where w = w{t), v = v{t), vi = v{ti). The next trick is important since it allows us to kill the "bad 
terms" [w{t) — v{ti)] ■ V in the operator A(t) — A\. 

Definition 7.1 Let us change the variables {y,t) >-^ {yi,t) = {y + di{t),t) where 

di{t) := {w{s) - v{ti))ds, 0<t<ti. (7.2) 

Jti 

Next define 

Zi{t) = {ei(yi,t),ai{yi,t),r{t),TT{t)) := {e{y,t),a{y,t),r{t),TT{t)) 

= {e{yi-di{t),t),aiyi-diit),t),r{t),7r{t)). (7.3) 

Then we obtain the final form of the "frozen equation" for the transversal dynamics 

Zi{t) = AiZi{t) + Bi{t)Zi{t) + Ni{t), 0<t<ti, (7.4) 

where Ni{t) = N{t) is expressed in terms of y = yi — di{t), and Bi{t) = 

/ Us{p{By^t) - By^t,^){p, •)) U,{-Vp{v{t) - v{ti))) -UMByit) - S.(to)-) \ 



-(-6^,(4) - i?^(tj) )(/?,•) By(^t)-By(^t,) 

\ {p, (vit) - viti)) ■ V-) -{-Vp, (f (t)V^(t) - viti)VAy(^t,^)) / 
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At the end of this section, we wih derive appropriate bounds for the "remainder terms" Bi{t)Zi{t) and 
Ni{t) in (j7.4p . First, note that we have by Lemma 16.21 



|i?.(t) - i?.(tol < I / Hs) ■ VvB,^s)ds\ < C / \\Z{s)\\tpds. (7.5) 

Jti Jt 

The similar estimates hold also for all the rest terms of Bi{t). Let us recall the following well-known 
inequality: for any a G IR 

(l + |2/ + x|r<(l + |2/|r(l + |x|)H, x,yGlR3. (7.6) 



Lemma 7.2 For (e, a, r, vr) € £^a wii/i any a € IR the following estimate holds: 

||(e(yi-di),a(yi-di),r,7r)|U<||(e,a,r,7r)|U(l + |(ii|)H , di G IR^. (7.7) 

For proof see |20| . 

Corollary 7.3 The following hound holds 

||iVi(t)||/3 < \\Zim%{l + \d^{t)\f^ , < t < ti. (7.8) 

Indeed, applying the previous lemma twice, once for (3 and once for — /?, we obtain from (j6.22p that 

iiivi(t)ii;3 < (1 + w{t)\nN{t, z{t))\\p < (1 + \d,{t)\nz\tp < (1 + \d,{t)\fp\\z,{t)\tp . 



Corollary 7.4 The following hound holds 



\\Bi{t)Z^{t)\\p<C\\Z^{t)\\.p j {l + \di{T)\f^\\Zi{T)f_pdT , 0<t<ti. (7.9) 

For the proof we apply Lemma 17.21 to (j7.5p and to the bounds for all the rest terms of Bi{t). 

8 Integral Inequality 

The equation (|7.4p can be written in the integral form: 

Zi(i) = e^i*Zi(0) + / e^i(*-')[SiZi(s) + Ni{s)]ds, 0<t<ti. (8.1) 

Jo 

We apply the symplectic orthogonal projection Pi := Pv{ti) to both sides, and get 

PiZi{t) = e^''PiZi{0)+ f e^i(*-^)Pi[5iZi(s)+iVi(s)]ds. 

Jo 

We have used here that Pi commutes with the group e^* since the space Zi := Pi£ is invariant with 
respect to e^* by Proposition 16.81 iil (cf. also Remark 16. 7p . Applying (|6.28p we obtain that 

||Pi^i(t)||-/3 < ^^-^||PiZi(0)||^ + cj^ __l__\\p^[B^Z^{s) + N^{s)]yds. (8.2) 
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The operator Pi = I — IIi is continuous in £^ by (j6.25p . Hence, from (j8.2p and (|7.8p . (j7.9p . we obtain 
that 



\\PiZi{i)\\^p 
+C(di{ti)) 



* 1 



ii + \t-s\r 



Zi(.)||_^ / ||Zi(r)rVT+||Zi (5)111^ 



ds, 0<t <ti, 



(8.3) 



where di{s) := supQ<«^ |(ii(t)|, < s < ti. Since \\Zi 
rewrite (18.31) as 



< C((ii)||Z(t)||_^ by Lemma [721 we can 



\PiZi 



+C{di 



Cidi] 

1- 
1 



-^^(TWll^^^^ll^ 



l^(s)|L 



(l + |t-s|)/3 



ti 



|Z(r)||V^+||Z(s)f„^ 



ds, 0<t <ti, 



•4) 



We can replace here the constants C{di) by C if di is bounded for ti > 0. In order to do this replacement, 
we reduce the exit time 



Definition 8.1 t'^ is the exit time 

t'^ = sup{t G [0,t*) : di{s) < 1, < s < t}. 

Now ([8^ implies that for ti < t'^ 

C 



(8.5) 



|Pi^i(t)||-/3<^^^P^||^(0)||^ 



+C 



1 



/o {l + \t-s\Y 
9 Symplectic Orthogonality 



Z{s)\\^p j'^' \\Z{T)\tpdT + \\Z{s) 



ds, 0<t< ti. 



(8.6) 



Finally, we are going to change PiZi{t) by Z{t) in the left hand side of (|8.6p . We will prove that it is 
possible using again that dp <^1 \n (12. 9p . For the justification we reduce the exit time once more. First, 
we introduce the "majorant" 



m{t):= sup(l + s)^||Z(s)||_^ 
se[o,i] 



t G [0,t=,). 



(9.1) 



Let us denote by e a fixed positive number which we will specify below. 
Definition 9.1 t" is the exit time 

t'l = sup{t E [0, t'J : m(s) < e, < s < t}. 



(9.2) 



The following important bound (j9.3p allows us to change the norm of PiZi(t) in the left hand side 
of ()8.6p by the norm of Z{t). 



Lemma 9.2 For sufficiently sm,all e > 0, we have for ti < t" 

||Z(t)||_^<C||PiZi(t)||_^, 0<t<ii, 

where C depends only on p and v. 



(9.3) 
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Proof The proof is based on the symplectic orthogonahty (j6.27p . i.e. 

n„(i)Z(t) = o, te[o,ti], 



(9.4) 



and on the fact that all the spaces Z{t) := Py{t)£ are "almost parallel" for all t. 

Namely, we first note that ||Z(t)||_^ < C||2'i(t)||_^ by Lemma[721 since \di{t)\ < 1 for t < ti < t'i < 
t'^. Therefore, it suffices to prove that 



iZi 



.^ < 2||PiZi 



0<t<ti. 



This estimate will follow from 



since PiZi{t) = Zi{t) - n„(j^)Zi(t). To prove ^M), we write ([93]) as 



n 



v{i),l 



Zi(t)=0, tG[0,ti], 



(9.5) 



(9.6) 



(9.7) 



where n^(t) iZi(t) is n^(t)Z(t) expressed in terms of the variable yi = y + di{t). Hence, ()9.6p follows 
from ()9.7p if the difference Yi^i^n) ~ ^v{t),i is small uniformly in t, i.e. 



|n„(t^)-n,(t),i|| <i/2, o<t<ti. 



(9.8) 



It remains to justify (|9.8p for a sufficiently small e > 0. 

In order to prove the bound (j9.8p . we will need the formula ()6.25p and the following relation which 
follows from (j6.25p : 

~ (9.9) 



n 



v(t),l 



Zi{t) =J2^jMt))r,,i{v{t))n{n,i{v{t)),Zi{t)), 



where Tj^i{v{t)) are the vectors Tj{v{t)) expressed in the variables yi. In detail (cf. ()3.4p ). 

j = 1,2,3, 



Tj^i{v) := {-djE^{yi-di{t)),-djA^{yi-di{t)),ej,0), 
Tj+3^i{v) := {dy^Ey{yi- di{t)),d^^A^{yi-di{t)),0,dy^Py), 



(9.10) 



where v = v{t). Since |(ii(t)| < 1 and Vtj are smooth and sufficiently fast decaying at infinity functions. 
Lemma 17.21 implies that 



for ah j = 1,2, 



and therefore 



Similarly, 



\h,Mt)) - rj{v{t))y < C\di{t)f, 0<t<ti 
, 6. Furthermore, 



(9.11) 



rMt)) - rMh)) 



ti 



vis) ■ VvTj{v{s))ds, 



\Tj{v{t)) - Tj{v{tl))y <C r \v{s)\ds, 0<t<ti. 



\n^i{v{t)) - Uji{v{ti))\ = \ r v{s) ■ V^Uji{v{s))ds\ <C r \v{s)\ds, 0<t<ti, 

Jt Jt 



(9.12) 



(9.13) 



since \V v^ji{v{s))\ is uniformly bounded by (|6.9p . Hence, the bounds (|9.8p will follow from (j6.25p . 
and (J9.1ip - (j9.13p if we establish that |di(t)| and the integral in the right hand side of ()9.12p can be made 
as small as we please by choosing e > sufficiently small. 
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To estimate di{t), we note that 

w{s) — v{ti) = w{s) — v{s) + v{s) — v{ti) = c{s) + j v{T)dT 
by (|6.10p . Hence, (j7.2p . Lemma [612] and the definition (|9.ip imply that 
\di{t)\ = 



(9.14) 



{w{s) - v{ti))ds\ < j |c(s)| + / \v{T)\dT]ds 



h 



since ti < t'l. Similarly, 



t K{l + sYf 



+ 



dT 



\v{s)\ds < Cm^{ti] 



(1 + t)2/3 



ds 



ds < Cm'^iti) < Ce^ 0<t<ti (9.15) 



(l + s)2/3 



< Ce^ 0<t<ti. 



The proof is completed. 



(9.16) 

D 



10 Decay of Transversal Component 



Here we prove Proposition 16.41 

Step i) We fix e > and t'l = t"{e) for which Lemma [9^2] holds . Then the bound of type (|8.6p holds with 

||PiZi(t)||_^ in the left hand side replaced by ||Z(t)||_/} : 



\Z 



+c 



-/3 



< 



c 



[l+tf^ 

1 



Z{s)\\.f,l^'\\Z{T)\\l^dT + \\Z{s) 



ds, 0<t<ti 



/o {i + \t-s\r 

for ti < t^. This implies an integral inequality for the majorant 

m{t) := sup {l + sf\\Z{s)\\-fs . 

sG[0,t] 

Namely, multiplying both sides of (jlO.ip by (1 + 1)^, and taking the supremum in t G [0,ti], we get 

m(s) Z"*^ m?{T)dT m?{s) 



(10,1) 



for fi < t'l- Taking into account that m{t) is a monotone increasing function, we get 



ds 



m(ti) < C\\Z{0)y + C[m\ti) + m\ti)]I{ti 



ti < tl 



(10.2) 



where 



fi (l + t)'^ 
'"'* = , S,io (1 + 1*-. I)" 



*i dr 1 

+ 



(1 + S)/3A (1+t)2/3 (1 + s)2/3 



ds < I < oo, ti > 0. 



Therefore, ()10.2p becomes 



m(ti) < C||Z(0)||;3 + C7I[m3(ti) + m\ti)], ti < t'l. 



(10.3) 
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This inequality implies that m{ti) is bounded for ti < t'l, and moreover, 

m(ti)<Ci||Z(0)||^, h<t:, (10.4) 

since m{0) = \\Z{0)\\^ is sufficiently small by (j3.7p . 

Step ii) The constant Ci in the estimate (|10.4p does not depend on t^, t'^ and t" by Lemma 19.21 We 
choose d/3 in (j2.9p so small that ||Z(0)||/3 < e/(2Ci). It is possible due to (j3.7p . Then the estimate (jl0.4p 
implies that t" = t'^ and therefore (|10.4p holds for all ti < t'^. Then the bound (I9.15P holds for all t <t'^. 
We choose e so small that the right hand side in (j9.15p does not exceed one. Then t'^ = i*. Therefore, 
(jl0.4p holds for all ti < t*, hence (16. 4p also holds if ||Z(0)||/3 is sufficiently small. Finally, this implies 
that t=K = oo, hence also t" = t'^ = oo and (jl0.4p holds for all ti > if dp is small enough. D 

11 Soliton Asymptotics 

Here we prove our main Theorem 12.21 under the assumption that the decay (j6.19p holds. First we will 

prove the asymptotics (|2.10p for the vector components, and afterwards the asymptotics (12. lip for the 

fields. 

Asymptotics for the vector components From (14. 2p we have q = b + r, and from (I6.2ip . (j6.22p . 

(frisj) it follows that r = -B^^t){p,a) + 5i,(t)vr + OdlZJI^^^). Thus, 

q = b + r = v{t) + c{t) - B,^t){p, «) + 5.(t)vr + 0(||Z||%). (11.1) 

Recall that (3 = 1 + 6, 0<d < 1/2. The equation (fOT]) and the estimates (f6T2]) . (f6T9]) imply that 



Ci{p,v, di3) 



m\ + \Ht)\< ,r'\',zi , t>o. (11.2) 



Therefore, c{t) = c+ + 0{t-^^+'^^'>) and v{t) = v+ + 0{t-^^+^^'>), t -^ cx). Since ||a||_i,/3 and |7r| decay like 
(1 + 1)-^, the estimate (|^T9]l . and (HOD, (flTT]) imply that 

q{t) = v+ + 0{t~^). (11.3) 

Similarly, 

b{t) = c{t)+ v{s)ds = v+t + a+ + 0{t-^^), (11.4) 

hence the second part of (I2.10p follows: 

q{t) = b{t) + r(i) = v+t + a+ + ©(r^-^), (11.5) 

since r{t) = 0{t~f^) by (I6T9D . 

Asymptotics for the fields We apply the approach developed in [18], see also [HI [171 (El [24]. For 

the field part of the solution, F{t) = {E{x,t),A{x,t)) let us define the accompanying soliton field as 
-^v(t)(*) = {Ev{t)ix ~ Qi't))^ ^v(t)ix — Q{t))), where we define now v(t) = q{t), cf. (|ll.ip . Then for the 
difference Z{t) = F{t) — Fv(t)(t) we obtain easily the equation [21], Eq. (2.5), 

Z(t) = AZ{t) - V • VvF,(i) (t), A{E, A) = i-AA, -E). 

Then 

Z{t) = W\t)Z{{)) - f W^{t - s)[v(s) • VvF,(,)(s)]ds. (11.6) 



'0 

To obtain the asymptotics (|2.1ip it suffices to prove that Z{t) = W^{t)"^^ + r^{t) with some ^+ G T 
and ||?'+(i)||jc- = 0{t^ ). This is equivalent to 

W\-t)Z{t) = ^+ + r'^{t), (11.7) 
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where ||r^(t)||jp = 0{t ) since W'^{t) is a unitary group in the Sobolev space J- by the energy conserva- 
tion for the free wave equation. Finally, (jll.Tp holds since (jll.6p implies that 

W'^i-t)Z{t) = Z{0)+ [ W\-s)R{s)ds, i?(s) = v(s)-VvF^(,)(s), (11.8) 

Jo 

where the integral in the right hand side of (Ill.Sp converges in the Hilbert space J-' with the rate 
0{t-^). The latter holds since \\W^{-s)R{s)\\jr = ©(s"^) by the unitarity of VF°(-s) and the decay 
rate ||i?(s)||:r = 0{s~^) which follows from the asymptotics for the vector components. More precisely, 
differentiating the equation (jl.3p in t and using the asymptotics (|11.3p . (|6.19p we obtain an estimate for 
v(t) = q{t) providing the mentioned decay rate of R{s). D 

12 Decay for the Linearized Dynamics 

In remaining sections we prove Proposition l6.8l in order to complete the proof of the main result (Theorem 
12. 2p . Here we discuss our general strategy of the proof of the Proposition. We apply the Fourier-Laplace 
transform 

f'OO 

X{X)= e-^^X{t)dt, ReA>0 (12.1) 

Jo 

to ()6.23p . According to Proposition 16.81 we expect that the solution X{t) is bounded in the norm || • ||_^. 
Then the integral (jl2.ip converges and is analytic for Re A > 0, and 

ll^(A)||-/3 < j^, ReA>0. (12.2) 

Let us derive an equation for X{X) which is equivalent to the Cauchy problem for ()6.23p with the initial 
condition X(0) = Xq G £~i3- We will write A and v instead of Ai and vi in all remaining part of the 
paper. Applying the Fourier-Laplace transform to ()6.23p . we get that 

XX{\) = AX{X) + Xo, ReA>0. (12.3) 

Let us stress that (|12.3p is equivalent to the Cauchy problem for the functions X(t) G C[0,oo;£-p). 
Hence the solution X{t) is given by 

X{X) = -{A- xy^Xo, ReA>0 (12.4) 

if the resolvent R{X) = {A — X)~^ exists for Re A > 0. 

Let us comment on our following strategy in proving the decay (|6.19p . First, we will construct the 
resolvent R{X) for ReA > and prove that it is a continuous operator in <S_^. Then X{X) G E-p and 
is an analytic function for Re A > 0. Second, we have to justify that there exist a (unique) function 
X{t) G C[Q,oo;£^p) satisfying (fT2l^ . 

The analyticity of X{X) and Paley- Wiener arguments (see [23J) should provide the existence of a i?_^ 
- valued distribution X{t), t G IR, with a support in [0,oo). Formally, 

X{t) = — I e''^^X{iuj + 0)duj, t G IR. (12.5) 

However, to check the continuity of X{t) for i > 0, we need additionally a bound for X{ioj + 0) at large 
|ci;|. Finally, for the time decay of X{t), we need an additional information on the smoothness and decay 
of X{iijj + 0). More precisely, we should prove that the function X{iijj + 0) 

i) is smooth for liJ 7^ 0, 

ii) decays in a certain sense as |a;| — > 00. 
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iii) is analytic at to = if Xq £ Zy := P^S and Xq £ £p. 
Then the decay (j6.19p would follow from the Fourier-Laplace representation (jl2.5p . 

However, we will check the properties of type i)-iii) only for the last two components f(A) and 7r(A) 
of the vector X{X) = (e(A), a(A),f(A),7f(A)). The properties provide the decay (|6.19p for the vector 
components r{t) and 7r(t) of the solution X{t). 

At the end, we analyze the field components e(x, t) and a{x, t) without the Fourier-Laplace transform. 
Namely, we prove the decay (16.19P for the field components directly from the field equations of the system 
()6.23p . using the decay of the components r(t),7r(t) and a version of strong Huygens principle for the 
Maxwell equation. 



13 Constructing the Resolvent 



First let us make a change of variables in the equation ()6.23p which will simplify significant the structure 
of the resolvent. The equation (J6.23P reads 



e 
d 
r 

TT 



V ■ Ve — Aa -|- Iis{f " Vpv — pB^^ir — {p,a))), 

—e + V ■ Va, 

Bv{t^ - {p,a)), 

{p, V ■ Va) - (r • Vp, V ■ VA^). 



(13.1) 



Put (f = TT — {p,a). Then tt = if + {p,a). If we prove a decay of ip and a, then vr has the corresponding 
decay as well. Further 

(f = TT - {p,d) = TT - {p,-e + v ■ Va) = {p, e) - {r ■ Vp, v ■ VA^) 

by the last equation of (jl3.ip . Thus, the system (jl3.ip is equivalent to the following system, 

e = V -Ve — Aa + Ils{r -Vpv — pByip), 

a = —e + V ■ Va, 

r = B^Lp, 

if = {p,e) - {r-Vp,v-VA^). 

Denote by the same letter A the operator 



(13.2) 



A 



/e\ / v-Ve- Aa + Us{r ■ Vpv - pB^cp) \ 

a _ —e + V ■ Va 

r ' B^ip, 

\ipj \{p,e)-{r-Vp,vVA,) J 



(13.3) 



Below we construct and study the resolvent of A. 

To justify the representation (jl2.4p . we construct the resolvent as a bounded operator in (?__^ for 
Re A > 0. We shall write (e(y), a{y),r, if) instead of (e(y. A), a(y, A), f(A), i^(A)) to simplify the notations. 
Then (fT23|) reads 

/ e \ / eo \ 



(^-A) 



a 
r 



\^ J 



ao 
\ V J 
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This gives the system of equations 

V ■ Ve — Aa + Ils{r ■ V pv — pB^ip) — Xe = —cq 

—e + V ■ Va — Aa = — qq 

Bjjip - Xr = -ro 

{p,e) - {r-Vp,v-VA^) - Xip = -ipo 

Step i) Let us consider the first two equations. In the Fourier space they become 

—i{kv)e + k'^a — Ils{i{kr)pv + pB,uip) — Xe = — eo 

— e — i{kv)a — Xa = — oq 

From the last equation we have 

e = — (A + i{kv))a + gq. 

Substitute to the first equation of (jl3.5p and obtain 

a = — ((A + ikv)ao - cq + IT), 11 := pUs{i{kr)v + B^ip)), 



where 

It is easy to see that 

Finally, 



b = b{X) = k'^ + {X + ikvf. 

b{X) /O for ReA >0. 
fe^oo + (A + ikv)eo - (A + ikv)fl 

b 



(13.4) 



(13.5) 
(13.6) 

(13.7) 

(13.8) 
(13.9) 

(13.10) 



Let us proceed to the fourth equations of (|13.4p . The equation reads 

{p, e) - {r ■ Vp, V ■ VAy) - Xtp = -ipo. 

From now on we use the system of coordinates in x-space in which v = (|t;|,0, 0), hence vk 
(jl3.10p and a straightforward computation we obtain 



\v\ki. By 



where 



$ = $(A,p,ao,eo) : 



{p, e) = '^ + Cr + Fif, 

{k'^ao + (A + ikv)eo)p 



D 



dk. 



Remark 13.1 Note that 

^ = Ft-.x{W\t){eo,ao),p), (13.11) 

where Ft^x is Laplace transform in t, VF^(t) is the first component of the dynamical group W{t) defined 
below by flOl) . 
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Further, 

C = .H| o' c. Vc, = -/ ^^(^ + -^;'-')'^^ l-g)d., (13.12) 

C3 / ^ ^(^) 

— , ' dk, j = 2, 3, (13.13) 



F= I 0^' Ih V /,= /i^±#^M(S_i)rf,, , = 1,2,3; (13.14) 



vj, 



3 



Z)(A) A;2 



recall that v = Vl — v^ ■ Further, — (r • V/?, w • VA„) = Gr, where 

By the change of variables k2 i— > A:3 we obtain that C2 = C3, f2 = fs, 92 = 93- Let us denote c := C2 = C3, 
/ := /2 = /3, g '■= 92 = 93- Note that ci + C2 + C3 = and thus, ci = — C2 — C3 = —2c. Similarly, 
5fi = —2g, and the matrices C and G simplify to 

/ -2c \ / -25 \ 

C = i|w| co\,G = v^lo 9 \. (13.16) 

V c y V Off/ 

Finally, the fourth equation becomes 

iC + G)r + {F-XE)^ = -{ipo + <^). (13.17) 

We write the third equation of (jl3.4p and the equation (J13.17P together in the form 



where M{X) is the 6 x 6-matrix 



*^<''i^)= -'(■;+*) • '"■^'*> 



*^(^' = ( c'+a f- XE ) ■ <i"^) 



Assume for a moment that the matrix M(A) is invertible for Re A > (later we shall prove that this the 
case indeed). Then 

^;)=M-(A)(_J;^J, ReA>0. (13.20) 

Finally, formulas (jl3.7p . (|13.10p . and (|13.20p give the expression of the resolvent R{X) = (A — X)^^, 
Re A > 0, in Fourier representation. 

Step a) Let us now proceed to the x-representation. We invert the matrix of the system (113. Sh and 
obtain 

-(ivk + X) k"^ \ ,,. . ..0 ,2i-i f —iivk + X) —k"^ \ 

-1 -{ivk + X) J LV I ; t J \^ ^ _(^y/^ + A) J 

Taking the inverse Fourier transform, we obtain the corresponding fundamental solution 

Gx{y) = {^'''\~^ ^.^_;^)ffA(y), (13.21) 
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where gx{y) is the unique tempered fundamental solution of the determinant 

D = D{X) = -A + {-vV + Xf. (13.22) 

Thus, 

Note that the denominator does not vanish for Re A > 0. This implies 

Lemma 13.2 The operator G\ with the integral kernel G\{y — y') is continuous as an operator from 
H^{m^) e L2(_r3) fg if2(jj^3) ^ h^(^r3^ for ReA > 0. 

Thus, formulas ()13.5p and (jl3.2ip imply the convolution representation 

e = -(^vV - X)gx*eo- Agx*aQ + {vV - X)gx*'n 

a = -g\ * eo - (w • V - X)gx * eo + 5a * n 



(13.24) 



Step Hi) Let us compute gx{y) explicitly. First consider the case v = 0. The fundamental solution of the 
operator —A + A^ is 

9x{y) = -ryj- (13.25) 

Thus, in the case u = we have 

Gx{y-y') ' ^ ^ ^ ^ 



1 —A J 47r|y — y'\ 

For general v = (|f|,0,0) with |t;| < 1 the denominator in (jl3.23p . which is the Fourier symbol of D 
reads 

D{k) = k^ + {i\v\ki + Xf = (1 - v'^)kl + kl + kl + 2i\v\kiX + A^ 

i\v\X ,0 ,2 ,2 ^'^ 

-^f + kl + kl + -— 
Set 7 = (1 — t;^)~^' ^ and return to x-space: 



{I - v^)iki + j^f + k^ + kl + ^. (13.26) 



D = -4(Vi + 7'b|A)2 -Vl-Vl + (7A)2. (13.27) 

Define yi := 72/1 and Vi := d/dyi. Then 

D = -(Vi + -i\v\Xf - vi - V^ + (7A)2. (13.28) 

Thus, its fundamental solution is 

^-'■f\\y\-^\v\Xyi 
9x{y) = ^— pq , y := (7^1,^2,2/3), (13.29) 

where Re A > 0. Let us note that 

0<Re7|v|A <Re7A, (13.30) 

since \v\ < 1. 

Step iv) Let us state the result which we have got above. 
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Lemma 13.3 i) The operator D = D{X) is invertible in L^(]R^) for Re A > and its fundamental solu- 
tion \13.29]) decays exponentially. 

a) Formula \13.29\) implies that, for every fixed y, the Green function g\{y) admits an analytic continu- 
ation (in the variable X) to the entire complex plane. 

Lemma 13.4 The matrix function M (X) (M~^(X)) admits an analytic (meromorphic) continuation from 
the domain Re A > to the entire complex plane. 



Proof The analytic continuation of M{X) exists by Lemma 113.31 ii) and the expressions (J13.12p - (jl3.14p 
since the function p[x) is compactly supported by (11.121) . The inverse matrix is then meromorphic since 
it exists for large Re A: this follows from (p39]) since C(A), F(A) ^ as Re A ^ oo by (pl^ - (fT3Ti]) .D 



14 Analiticity for Re A > 

Proposition 14.1 The operator-valued function R{X) : £ ^ £ is analytic for Re A > 0. 

Proof It is sufficient to prove that the operator A — X : £ ^^ £ has a bounded inverse operator for 
Re A > 0. Let us recall, that A = A^^^ where \v\ < 1. 

Step i) Let us prove that Ker (A — X) = for Re A > 0. Indeed, let us assume that Xx = (e, a, r,ip) £ £ 
satisfies {A — X)Xx = 0, that is Xx is a solution to (|13.4p with cq = oq = and rQ = (po = 0. We have to 
prove that Xx = 0. 

First let us check that ipx = 0. Indeed, the trajectory X := Xxe^^ G C(IR, i5) is the solution to the 
equation X = AX that is (15.13P with w = v. Then TCv,v{X{t)) grows exponentially by (15.19p . since the 
matrix B^ is positive. This growth contradicts to the conservation of TCv,v Now Xrx = B^ipx = by the 
third equation of (|13.4p . hence r^ = since A 7^ 0. Finally, ba = 0, oa = by the equations (113. 24p with 
r = ta = 0. 



Step ii) One has 
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Thus, A- X = Aq + T, where 
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The operator T is finite-dimensional, and the operator A^ is bounded in £ by Lemma 113.31 Finally, 
^4 — A = Aq[I -ir Aq T), where A^ T is a compact operator. Since we know that Ker (/ + A^ T) = 0, 
the operator (/ + Aq T) is invertible by Fredholm theory. D 

Corollary 14.2 The matrix M(A) of 1113.1^) is invertible for ReA > 0. 

15 Regularity on the Imaginary Axis 

Next step should be an investigation of the limit values of the resolvent -R(A) at the imaginary axis 
A = iuj, a; G IR, that is necessary for proving the decay ()6.19p of the solution X{t) = (e(t), a(t), r(t), (p{t)). 
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However, we consider only the limit values of the matrix M~^(A), inverse to (|13.19p . This allows us to 
prove the decay of the finite dimensional components r(t), (p{t) by the inverse Fourier-Laplace transform. 
We will derive the decay of the field components e{t),a{t) later without the Fourier-Laplace transform. 
First, let us describe the continuous spectrum of the operator A = A^^^ on the imaginary axis. By 
definition, the continuous spectrum corresponds to uj £ IR, such that the resolvent R{iuj + 0) is not a 
bounded operator in £. By the formulas (jl3.24p . this is the case when the Green function g\{y — y') loses 
the exponential decay. This is equivalent to the condition that Re x = where >c = 7V— ti?. Thus, iio 
belongs to the continuous spectrum for all u; G IR. By Lemma 113.4^ the limit matrix 

M{ioj):=Miioj + 0)=( ''"^^ nN ^ ^^ nN ■ c. V ^elR, (15.1) 

^ ^ ^ ' \ C{iuj + 0) -h G F{uv + 0) - iujE y ' ^ \ ) 

exists, and its entries are continuous functions of u; E IR, smooth for \ijj\ > 0. Recall that the point A = 
belongs to the discrete spectrum of the operator A by Lemma 15.31 i). hence M{iuj + 0) (probably) is also 
not invertible at u; = 0. 

Proposition 15.1 Let p satisfy the Wiener condition U.13\) . and \v\ < 1. Then the limit matrix M{iuj + 
0) is invertible for uj ^ 0, uj £ IR. 

Proof Since v = (|w|,0,0), the matrix By is also diagonal: 

^3 
By := u{E - V (g) v) = { u |. (15.2) 

By (fT5Tl) and (fTXT2l) - (fTXTHI) . for a; G IR, 

detM(iu;)=det f "''^^ hn ^ w- n^ ■ c^ ^ = A^^ 

^ ' \ C{ioj + Q) + G F{iuj + 0) - iujE J 

where 

Di = -uj"^ - iojiy^fi{iuj + 0) - iy^i\v\ci{iuj + 0) - v'^v^gi, (15.3) 

D = —to — iujvf{iuj + 0) — vi\v\c{iuj + ^) — vv g. (15.4) 

The formula for the determinant is obvious since all of the matrices C, F, G, and By are diagonal. 
For |i^| > the invertibility of M{iuj) follows from (jl5.3p . (jl5.4p by the following lemma. 

Lemma 15.2 // il.l3\) holds and \uj\ > 0, then the imaginary parts of Di and D are positive: ImDi > 0, 
lmD>0. 

Proof Let cj > 0, the case liJ < is similar. Note that 

ImDi = -iy^{ujRefi{iuj + 0) + |7;|Re ci{iio + 0)). 

For e > we have 

;,(.„ + ,)= /(!^i±£±*HM*£,| _!),,, , = !,,, 3. (15.5) 

J D{iuj + e,k) k^ 

By Plemelj formula for G^'^-functions, cf. [30] (see Chapter VII, formula (58)) and ^Oj, 

Re/,X.. + 0)=./ <^±|li»£(|-lMS, (15.6) 

28 



where 



ri / 1 \'^\^ \2 ,1 t2 ^^t 

{k : {uki - ^-^-Y + k^ + kl = —} 



is the eUipsoid on which D{iuj, k) = 0. Similarly, 

Then 

u;Re/i(.a; + 0) + \v\Re ci(.a; + 0) =7r /" (^ + ki\v\?\m\\ ^ _ ^^^g_ 

Jt^ \VD{iLO,k)\ k^ 

By the Wiener condition p.lSp we obtain that the last integral is negative and hence IniDi > 0. 
Further, 

ImD = -u{ujRef{iuj + 0) + \v\Re c{iuj + 0)). 

Recall that / := /2 = /s and by (jl3.14p we obtain that 

J D{iuj + e,k) 2fc^ 

Then 

Jt^ \VD{iu;,k)\ ^ 2k^ ^ 

^r {u + kMmk)\\{u; + k,\v\f-kf 
Jt^ \VD{iu;,k)\ ^ 2fc2 

since /e| + A;| = [to + ki\v\Y — kf on T^^. Now recall that c := C2 = 03 and ci = —2c, by (|13.12p we obtain 

that 

^ W k,i^u. + e + ^kM)m)\\kJ _ 

27 D{iuj + e,k) k^ 

Then 

We take into account that k'^ = {to + A;i|t'|)'^ on T^^ and obtain wRe f{iio + 0) + |w|Re c{iio + 0) = 

vr f \p{k)\Hu^ + h\v\)Hu;^ + kUl-v^^^ 



2Jt^ \VD{iLO,k)\k^ 

This value is negative, again by the Wiener condition, hence Im D > 0. This completes the proofs of the 
lemma and the Proposition 15.1. D 

Corollary 15.3 Proposition \15.i\ implies that the matrix M~^{iuj) is smooth in w G IR \ 0. 

Remark 15.4 The proof of the Lemma [15. 2 1 is the unique point in the paper where the Wiener condition 
is indispensable. 
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16 Singular Spectral Points 

Let us recall that the formula (|13.20p expresses the Fourier-Laplace transforms f(A),(^(A). Hence, the 
components are given by the Fourier integral 

^^5 ] =— [e'^^'M-^iuj + 0)( , ""^^ ^, ^ duj (16.1) 

if it converges in the sense of distributions. The Corollarv 115.31 alone is not sufficient for the proof of the 
convergence and decay of the integral. Namely, we need an additional information about a regularity of 
the matrix M~^{iuj) at its singular point to = and some bounds at |u;| — > oo. We will analyze all the 
points separately. 

I. First we study the asymptotic behavior of M^^{X) at infinity. Let us recall that M~^{X) was originally 
defined for Re A > 0, but it admits a meromorphic continuation to the entire complex plane C (see Lemma 

MM- 

The following Lemma is a very particular case of a general fundamental theorem about the bound 
for the truncated resolvent on the continuous spectrum. The bound plays a crucial role in the study of 
the long-time asymptotics of general linear hyperbolic PDEs, |40j. 

Proposition 16.1 There exist a matrix Rq and a matrix-function Ri{uj), such that 



M-^(iuj) = — + Ri(u;), Iwl > 1, wGlR, 



where, for every k = 0, 1, 2..., 

|aii?i(w)| < -%, |w|>l, toeU. (16.2) 

The proof is similar to that of Proposition 16.2 in [20], see also [Ml Thm 3], [H the bound (A.2')], [221 
ThmS.l], [391 Thm 3]). 

II. Finally, we consider the point w = which is most singular. The point is an isolated pole of a 
finite degree by Lemma 113.41 hence the Laurent expansion holds, 

n 

M-i(iu;) = ^Lfca;"'''"^ + /i((j), |w| < eo, (16.3) 

A;=0 

where L^ are 6x6 complex matrices, Eq > 0, and h{uj) is an analytic matrix-valued function for complex 
CO with |lij| < eq. 

17 Time Decay of the Vector Components 

Here we prove the decay (|6.19|) for the components r{t) and (p{t). 
Lemma 17.1 Let Xq G Z^^/^. Then Q{t), P{t) are continuous and 

\r{t)\ + \^it)\<^0^, t>0. (17.4) 

Proof The expansions (jl6.2p and (jl6.3p imply the convergence of the Fourier integral (jl6.ip in the 
sense of distributions to a continuous function of t > 0. Let us prove (|17.4p . First let us note that the 
condition Xq G Zy^/^ implies that the whole trajectory X{t) lies in Z^^jj. This follows from the invariance 
of the space 2^t,,/3 under the generator Ay^y (cf. Remark 16. 7p . Note that for Xq not belonging to Zy^i^ 

30 



the components Q{t) and P{t) may contain non-decaying terms which correspond to the singular point 
u; = 0. Indeed, we know that the Unearized dynamics admits the secular solutions without decay, see 
(|6.24p . The formulas (j3.4p give the corresponding components rs(t) and ips{t) of the secular solutions, 



rsit) 



E^. 







t + 







8 P 



+ E^. 







(17.5) 



We will show that the symplectic orthogonality condition leads to (jl7.4p . Let us split the Fourier integral 
(jl6.ip into two terms using the partition of unity C,i{ijj) + C2('^) = 1, w G M: 



r{t) 



1 

2^ 



e^-*(CiH + C2H)M-^(iu; + 0) 



-((^0 + ^) 



duj 



= hit) + hit), 
where the functions Cfc(w) £ C°°(IR) are supported by 

suppCi C {lv eM: eo/2 <\uj\} 

suppC2 C {liJ e IR : |w| < eo} 

Then 

i) The function /i(t) decays like of (1 + |t|)~^ due to Proposition 116.11 
First let us prove the continuity. Consider the contribution of 



(17.6) 



(17.7) 



Ci 



Rn 



UJ 



-ro 

-^0 



'--Y.F.- + Y.F.-)Ro 

' U) UJ UJ 



-ro 

-^0 



Here jj- — V.P. ^ forms the compactly supported distribution which inverse Laplace transform is smooth. 
The inverse Laplace transform of V.P. — is sgn t which is continuous for t > (in the sense of disributions). 
The contribution of 

^''^'<-' ( :;: 

is continuous, since CiRiiuj) is summable by Proposition I16.11 Further, consider the contribution of 



^R ( ° 

a; °V^t^-(^'(i)(eo,«o),p) 



^-1 C 



Note that F~^i^Rq is continuous by the previous argument and has the decay t~°° . The function 
{W^it)ieQ, ao), p) is continuous and decays like t~^ by Lemma fl 8. 3 1 Thus, we come to the convolution of 
two continuous functions with decays t~°° and t~^. This convolution is clearly continuous. The similar 
argument is applied to the term 



Cii?i(cu) 





Ft^^{WHt)ieo,ao),p) 



The continuity is proved. 

Now let us prove the decay like t~^. We need to consider only the contribution of the term 

Cii-Ro + Riiuj))( ~ /wi.°v ^ \ 
OJ V Ft^Lu{W\t)ieo,aQ),p) 
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All the rest terms give the decay t °°: differentiate any times in uj and obtain a summable function of 
Lj. Consider 

Ci— -Ro 



u;''\Ft^^{WHt){eo,ao),p) 

The inverse Laplace transform gives the convolution of functions with decays t~°° and t~^. This convo- 
lution deacys like t~" . The last term can be considered similarly. 

ii) The function /2(t) decays like i~°° if (cq, ao,ro, (/^q) G 2^^^. 

We claim that the operator A\zp^ is invertible (recall that Z/^^^ is invariant w.r.t. A). First let us 
show that it has the trivial kernel. We have to resolve the system (jl3.4p with {eo,ao,ro,ipQ) = and 
check that (e, a,r,Lp) = if y = (e, a,r,ip + {a, p)) satisfies the symplectic orthogonality conditions. The 
vector part of the system, ()13.18p becomes 

hence ip = 0. The equations (|13.7p . (jlS.lOp give a = ipIls{kr)v/D, e = —ikvpIls{kr)v/D. For j = 1,2,3 
the symplectic orthogonality conditions ^(Xj'Tj) = are satisfied automatically and imply no restriction 
on r. The symplectic orthogonality conditions ri(y, Tj+3) = transform to Mr = 0, where M is a positive 
matrix (see the details in Appendix C), and thus, r = 0. This imply a = 0, e = and the triviality of 
the kernel is proved. 

Further, A acts in the two invariant subspaces: in Zp^^ and in its 6D symplectic orthogonal comple- 
ment; and the index of A is zero since it is an elliptic operator. Consider a homotopy which takes A to 
zero on this 6D subspace. Then the index of A\z0 „ remains zero and A\zi^ ^ has the trivial kernel. Then 
it has the trivial co-kernel and thus is invertible. 

Thus, R{X)\zp^ is analytic for A = ito, \oj\ < £q. Now, 

M-\iuj + Q)( "'■° , 
V -(V50 + ^) 

is the vector component of R[\)\zu ,„(-^o) a-i^d hence is analytic in A or lo. Then the l2{t) decays like t~°° . 

D 

18 Time Decay of Fields 

Here we construct the field components e{x , t) , a{x , t) of the solution X{t) and prove their decay corre- 
sponding to (I6.19p . Let us denote F{t) = (e(-, t), a{-,t)). We will construct the fields solving the first two 
equations of (j6.23p . where A is given by ()4.15p . These two equations have the form 

m = ( %^ ~%^F+(^ ^(*) ) , n(t) := UMt)Vpv - pBMt)). (18.1) 

By Lemma ll7.1l we know that r{t) and ip{t) are continuous and 

Ht)\ + \^{t)\<^0^, t>0. (18.2) 

Hence, the Proposition 16.81 is reduced now to the following 

Proposition 18.1 i) Let functions r{t),ip{t) € C[0,cxd;]R ), and Fq G J^. Then the equation 1118. 1\} 

admits a unique solution F{t) £ C[0, 00; J^) with the initial condition F{0) = Fq. 

ii) If Fq G J^p and the decay 1118.2]) holds, the corresponding fields also decay uniformly in v: 

\\m\\-P<^0^\\Foh, t>0, (18.3) 

for \v\ < V with any v £ (0, 1). 
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Proof The statement i) follows from the Duhamel representation 



Fit) = Wit)Fo + 



yit-s)^^^'^ ]ds 



t > 0, (18.4) 



where W{t) is the dynamical group of the modified wave equation 



The group W{t) can be expressed through the group Wo{t) of the wave equation 

6(i)=( _° "^ )$(t). (18.6) 

Namely, the problem (118. 6p corresponds to ()18.5p . when v = 0, and it is easy to see that 

[W{t)F{0)]{x) = [Wo{t)F{0)]{x + vt), X G IR^ t e M. (18.7) 

The identity (118. Th implies the energy conservation law for the group W{t). Namely, for 
(e(-,i),a(-,i)) = W{t)F{0) we have 

/ (|e(x,t) -u • Va(x,t)p + |Va(x,f)p)dx = const, i € M. 

In particular, this gives that 

\\W{t)Foy<C{v)\\Foy, tel. (18.8) 

This estimate and (jl8.4p imply the statement i). 

Let us proceed to the statement ii). Step i) 

Lemma 18.2 For any Fq £ Tfj the following decay holds, 

||W^o(t)i^o||-/3<^^[^||i^o||/3, t>0, (18.9) 

for the dynamical group Wo(t) corresponding to the wave equation il8.6\) . 

For proof see Appendix B. 

Step ii) 

Lemma 18.3 For v < 1 and Fq G Fp, the following decay holds, 

\\W{t)Fo\\.p < ^^\\Foy, t>0, (18.10) 

for the dynamical group W{t) corresponding to the modified wave equation il8.5\) with \v\ < v. 

Proof Let {e{x,t),a{x,t)) = [W{t)Fo]{x), {e{x,t),a{x,t)) = [Wo{t)Fo]{x). Since [W{t)Fo]{x) = [Wo{t)Fo] 
(x + vt), we have e{x, t) = e{x + vt, t), a{x, t) = d{x + vt, t). The norms ||a(-, t)||i^_^ and ||e(-, t)||_/3 are 
bounded by C(/3)||Fo||/3(l + |t|)-^. Let us start with e: 

/(I + \x\r^Vi^, t)dx ^ /(I + \x\r^^mx + vt, t)dx ^ /(I + I. - ..|)-^|ef (., t)dy < 

C{v) j{l + \y\)-''^\~e\\y,t)dy, 
33 



since |f| < f < 1. The last integral is bounded by C(/3)||Fo||o(l + \t\) ^^. The norm ||a||i^_/3 = 
11(1 + |x|)~'^Va||/,2 can be estimated similarly. 

Step in) Now the statement ii) follows from the lemma 118.31 and the Duhamel representation (jl8.4p . 
Indeed, W{t)Fo decays like t~^ by lemma [TS. 31 Further, 

t t 

W f W{t- s){j^^^^ ys\\.p< j \\W{t- s){j^^^^ y\^pds < 



} \ms)\\.pds ^ I ds 

J {i + {t-s)r- J {i + {t-s)ni+sr 



by Lemma 118.31 and (|18.2p . The last integral decays like (1 + 1)~^ by a well-known result on decay of a 
convolution. D 

A Appendix: Computing Symplectic Form 

Here we compute the matrix elements Q{tj,ti) of the matrix and prove that the matrix is non- 
degenerate. For j, / = 1, 2, 3 it follows from (|3.4p and (|3.2p that 

n{Tj,Ti) = {djE,, diA,) - {djA,,diE,), (a.i) 

n{Tj+3, T1+3) = {dy^E^, ^y^A^) - {d^^A^, dy^E^), (A. 2) 

and 

il(Tj, r^+s) = -{djEy.d^^A^) + {djA^.d^^Ey) + Cj • dy^P^. (A.3) 

In Fourier representation the solitons read 

E.m = ^ (M, _ „) , AAk) = ^ (^^ - ..) . (A,4) 

, . . f \p\^dk f \p\^dk,, ,, ,. , 

Py=Py + {Ay, p)=py + v j ^^ - / yjL_^kv)k, (A.5) 

where D := k^ — (kv)'^; D is nonnegative and even in k. Differentiating in v we obtain for j = 1, 2, 3: 

a.,4 = g(^^- '=-"'^;'^-"'^ -(Me,), (A.6) 

1 f\p\'dk f\pl{kv)kidk f mi^jHkvlJkdk 
Then for j, / = 1, 2, 3 we get from (jA.ip by the Parseval identity, 

{djEy,diAy) = -i ' '^^/'^' ( ^k -vj dk = 0, 

since the integrand function is odd in k. Similarly, {djAy, diEy) = and thus 

n{Tj,Ti) = 0. (A.9) 
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Further, by I^KM, 

{OvjK^dy^A^} -I j — I — + 

2kjki{kv){k'^ + {kvf) 2kjki{kv){k'^ + {kvf)v^ kjik"^ + {kvf)vi 



D2 1)2 JJ 

^^-^-^^- '-('"^'^' "MMS ikv)6,, + '^'--- ]dk = 0, 



k^kAkvMk"^ + (kvff 2ki{kvfv. , , k.ki(kv)(k^ + (kv'^^ 
1^2 D ^^")^^' + k^ 



since the integrand function is odd in k. Note that the integral converges by the neutrahty condition 
(I1.14p . Similarly, {dvjAy,dviEy)=0 and thus, 

fi(r,+3,ri+3)=0. (A.IO) 

Now let us compute Q{tj,ti^3). First, 

/f,^.^A\ fk,{kv)\p\^dk/ 2k{kvy ki{kv){k'' + {kvf) 

-{Ojt!jv,OviAy) = / —2 Vi + 



D^ \' D k^D 

Second, 

m A ^ i7\ fblP^ffi. ^ , k{k^ + {kv)^y 2ki{kvf \ 



And third, 



1 , i^dk f mkv)kidk f \p\\e + {kvf)k,kdk 

Cj ■ o^ r-y — Cj ■ n^ ei -\- Oji -\- ^ „2 "^i 



D ^' J Z)2 "^ J k^D^ 

We use the rotational invariance of p and apply the orthogonal change of variables k = Vk', v = Vv\ 
where v' = (|v|,0,0) and V = \\Vji\\ is an orthogonal matrix with Vji = Vj/\v\. Then 

y{kv)kidk _ f \y{k'v'){Vk')idk' _ f \y{k[\v\){Vk')idk 



D^ J (A;'2 - (A:V)2)2 J ^k'^ _ kfy2)2 ■ 

In the expression /c'^^juKV^iA;^ + Vi2k2 + Visk'^) we keep only the term /c'l^lwlV^i = k'lVi, since all the rest 
integrals equal zero. Finally, the integral equals 



/ 



k?\ydk' 
vi 



JJ/2 

(here D' = k'"^ — k'^v^ in the new coordinates) and 

\p?'^k^ , ^ [ k'^\p?dk' _^ [\ykjkidk [ \p\'^{kvfkjkidk 



.,1 r\prdk^ fk['\p\'dk' f 

e, ■ d,,P, = e, ■ B-'ei + J ^-^^ji + 2 j ' '^',, v,vi - J ■ 



D^ .1 k^D^ 

In total, 

^, , 1 /■ l^dk . f k'?\p\'^dk' 

2 f \p\^k'^kjkidk f \p\'^{kv)^kjkidk 



L>3 J k'^D^ 

2^ f \p?{kv)%kidk f \p\^k,kldk f lyjkvfkjkidk 

-•^i^ ^V d'^ j d'^ j k^^ ■ 

(A.ll) 
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The matrix elements of the first hne are those of a diagonal positive and nonnegative matrices. In the 

second line we have 

9p(i;^fc^ - {kvY)kjkidk 

Since v'^k'^ — {kv)^ > 0, this is the matrix element kjki of the nonnegative matrix k k integrated with 
a nonnegative weight, and thus the result is the matrix element of a nonnegative matrix. 

Observe that in the third line the matrix element kjki is integrated every time with a non-positive 
weight. We apply the described change of variables to all of these integrals. Then we may omit the prime 
in all of them. We keep only the terms that lead to nonzero integrals and represent the most essential 
part of the change as the following chain of transforms: 

kjki ^ {Vjiki + Vj2k2 + Vj3k3){Viiki + Vi2k2 + Vi^ks) ^ VjiViikf + Vj2Vi2k1 + Vj^Vi^kl 

since all the rest terms lead to zero integrals. Since the integrals with A;| and A;| are equal by the 
corresponding change of variables, 

VjiVnkj + Vj2Vi2kj + VjsViskI ^ VjiVnkj + {VJ2V12 + FjsV/s)^:^ ^ 

VjiVnkl + {Vj2Vi2 + Vj3Vi3 + VjiVn)kl - VjiVnkj = VjiVnikf - k^) 
by the orthogonality of the matrix V. Further, 

V]iVii{Ki - K2) - —{f^i - 1^2) ^ ~y'^^ 2 ' ~ ~iJ2~v% 2 ' ~ ~iP 2 ■ 

Then 

\p\'^{kvfkjkidk VjVi f \p\'^kf{3kf - k'^)dk 



L»3 2 J D^ 

\p\'^kjkidk VjVi f |/9p(3A;f — k'^)dk 



Z)2 2^2 J D^ 

d\^{kvfkjkidk VjVi f \p\^kl{3kl - k^)dk 



BD"^ 2 J k'^D'^ 

Finally, in the formula (|A.lip for i7(rj,ri_|_3) we obtain VjVi multiplied by the following expression: 

9 3 1 3(l-t;^) 9(1 -t;^) 3 



where 



i,= /M!M, j2= /"^^, /3= /^S^, h= /Si, h- ^'^''^' 



Z?2 ' ^ / L)2 ' ^ / 1)3 ' * / D3 ' / 1^2JJ2- 



All of the integrals Ii to Is can be computed in terms of |u| (in further we write v instead, for simplic- 
ity of notation) and R := Jq °° \p\'^dr by the spherical change of variables. After the straightforward 
computation we get 

irR f I 1 ^ l-v\ ^ vri?/ 1+f 1 

A = -3- ( :; - ^— + log :;— ) , I2 = ^ ( log :; + 



V^ \1 — V 1 + V 1 + V J ^V 1 — "U 1 — V 1 + V 

irR / l-v 1 1 1 1 

l3 = -r^{ log —— + — + 



4f3V 1 + w 1 + v l-v (l-t;)2 (l + v)'^ 

ttR /^^ 1 + v 5 5 1 1 

7-? 3 log -— + -— - -— + 



4v^ \ ^l-v l + v l-v (l-t;)2 (l + vY' 
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inR ttR f ^ 1-v 1 1 



V^ v^ \ 1 + v 1 — v 1 + f^ 

It remains to make Taylor expansion in v and observe that the expression ()A.12p involves only the even 
powers of v: f ~^, v^, f^, v^, ... ; the coefficient at v~'^ is zero; the coefficients at all the rest powers of v 
are positive. As the result, the matrix 

fi+(7;) = ||f](T,-,rz+3)IIU«=i,2,3 (A.13) 

is positive definite and hence non-degenerate. Finally, the matrix 

ll^(r„rO||U,=i,...,6=( V(.) r^'O ^^-^^^ 

is also non-degenerate. 

B Appendix: Decay of the Group of Free Wave Equation in Weighted 
Sobolev Spaces 

We prove Lemma 118.21 for the wave equation written in the standard form 

n, .[ID m. nf> . ( ;(') ) . ,b., 

For simplicity of notations let us denote the corresponding group W{t) instead of W^{t). Denote by 

W{x — y, t) the (distribution) integral matrix kernel of the operator W{t). The Huygens principle for the 
group W{t) reads 

W{x -y,t) = 0, \x-y\> t. (B.2) 

The energy conservation law for the group W{t) has the form 

/ (|e(x,t)p + |Va(x,t)p) dx = const, t e IR. 

In particular, this gives that 

\\F{t)\\:F<C\\Fo\\:F, tGM. (B.3) 

The matrix kernel W{z, t) can be written explicitly since the solution to (IB.ip has the form 

a{;t) = [—R{t)*ao + R{t)*eo], e(-, t) = d(-,i). (B.4) 

Here R{t) = R{-,t), and 

Step ii) It remains to estimate ||VF(i)Foll-/3 ^or large t > 0. Let us fix an arbitrary t > 1, and split the 
initial function Fq in two terms, Fq = Fq^ + Fq^ such that 

\\n,th + \\Flth<C\\F4p, t>l, (B.6) 

where C does not depend on t, and 

F^,i(x)=0, \x\>t + l, (B.7) 

Fl^{x) = 0, \x\ < t. (B.8) 
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Now the estimate for W{t)F^^ follows by (jES]), (iBlSll and (IRB]) : 

< Ci||F^^J^(l + |t|)-^<C2||Fo||;3(l + |t|)-^, t>l. (B.9) 

iStep in) It remains to estimate W{t)FQ^. By ()B.4p . (jB.SP The support of e(-,f),a(-,f) is in the area 
{\x\ > t/2}. Then 

11(1 + \x\)-^e{t)\\L2 < C(l + t)-^\\emL^ < C\\Foh2 < C\\Foy. 
The norm ||(1 + |x|)~^Va(t)||j;^2 cand be estimated similarly. 



C Appendix: Computations for Invertibility of ^41^^ 

n{Y, Tj) = -{e, djA^) + (a, djE^) - {a, p) • ej 



■'I3,v 

For j = 1,2,3 we have 



dk e ^^^2^_±lQk -v)+ I dka --^H--^/- (-fc _ ^) _ („^ ^) . e^. 



-ikj{—p),kv^ , /" ,, ^ —ikd{kv)p,kv 
Since e±/c and a_Lfc, the expression simplifies to 



/dfc e ^^t; — I dk a — v — I dkdp ■ e,. 
D j D J ^ ' 



—i I dk e 

We substitute 

(kv)(kr)p kv ^ . ^ i(kr)p kv 

and obtain that the expression equals zero, since every integrand function is odd in k. 
Further, we take the relations e_Lfc, o_Lfc into account and compute 

n{Y, Tj+s) = (e, d^^A^) - {a, d^^Ey) + r ■ d^Ty = 

dk e j-i ^: ' v + ej] -I dk a j-i -^ ^-^ — ^—v + {kv)ej j + r • dy^P^. 

We substitute the expressions for e and a and obtain: 

\p\'^dk 



n{Y,T,+,) = r.B-'e, + l 



D ' 

2 / ^p\^{kvfdk.^ ,, f\p\^kv)^dk.^ ., 



^3 {kr)k, - V / ^3 {kr)kj - 3 / ^^^^ ' 



'^' {kr)k, + 2 / '^' '^' {kr)v, + 2 / \!l^^—{rv)kj. 



j^2 ^ ' J ' J JJ2 '^ ' J ' J D^ 

In the first line there is the j'-th component of Mir, where Mi is a positive matrix. The second line 
equals 

pWkvy{k'^v^-{kv)^)dk 



Note that k'^v'^ — (kv)'^ > 0. Thus, one applies the nonnegative matrix k(>i)k integrated with a nonnegative 
weight to r and takes the j'-th component. To all of the rest integrals we apply the change of variables 
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described in Appendix A. We make this change, omit the prime, and keep only the terms that lead to 
nonzero integrals. Thus, 

3 3 3 3 

{kr)k, ^ {Vk ■ r){Vk)j = [(^ Vuk)n + (J^ ^2(^0^2 + {J2 ^s/^O^s] ^ Vjik ^ 

1=1 1=1 1=1 1=1 

3 3 3 

(E ViiVjikf)n + (Y^ V2iV,ikf)r2 + (^ V3iV,ikf)r; ^ 
1=1 1=1 1=1 

{ViiVjikf + {VuVj2 + I^i3^i3 + VnVji)kl - VuVjikDn + ■■• = 

VlVjfj2 ,2x , V2Vj , 2 ,2n , ^3^'j.,2 , 2x ^j-k^ 

-^ih - hVi + -^{ki - k2)r2 + -^{ki - fcajrs = ^ — ((^ ^ ^)'^)j- 



Further, 
and finally 
Then 



yZ yZ yZ yZ 



1,2 I 1.2 i.2 I 1.2 I f,2 _ t.2 ot.2 _ ^,2 

„2 u2 . , ,^2 '^2 T^ '*'3 _ ,2 '^2 T^ '^3 T^ '^1 '^1 _ "^'^1 ''' 



"-1 "'2 '*'! o — '''1 



{kr)kj h^ — {{v (E) v)r)j. 

\pWkv)^dk 1 f \p\'^kl{3kf - k'^)dk 



^ {kr)kj ^ - I Pig {{v v)r)j, 



Z?3 ^■"■^■"^ 2 7 Z)3 

^^^k'^dk 1 f\p\'^k'^{3kj-k'^)dk 



D^ ' ' ' 2 J Z)3 

/5|2(A:t;)2dA: 3 f |p|2A;2(3fc2 _ A;2)(iA; 



(A:r)fcj ^ - / -r— {{v®v)r)j, 



1^2 JJ2 y^ ' :> 2 J k^D^ 

\p\'^dk 1 f Ip^Skl - k'^)dk 



D2 (^O^j' ^ ^ y ^2 ((^ ® ^)^)j- 

Further, by the similar argument 

ki ki 

{kr)vj >-^ TT^^^ '^ ^)'')i ^^d {rv)kj >-^ 77 ((^ ® ^)'") 

Then 

|/5|2(A;i;)dA:,, , f \p\'^kUk ,, 

\2(U.,\^U (■ \?.\2u2. 



]■ 



\p\\kv)dk^ ,, ^ [ Wkidk^^ 



In total we obtain {{v ® v)r)j multiplied by the following expression: 

11 3 ,, 1 , 3, 9, 1 

2 2t;2 ^ ^ ^ 2^2 2 3 T 2 4 2^2 

where the integrals /i to I5 are introduced in Appendix A, and 



— - ^^2)^1 + TT^2 - 2/3 + 7^4 - -/5 + 7^6, (C.l) 



pfc^dA; ttR / l + v 6v 1 

^ 31og- + - ^ + 



L>3 4„ \^ "'l-t; l-i;2 (l-w)2 (1+?;)^ 

It remains, similar to Appendix A, to make Taylor expansion in v and observe that the expression (jC.ip 
involves only the even powers of v: v~'^, v^, v"^, v'^, ... ; the coefficient at u~2 ig zero; the coefficients at 
all the rest powers of v are positive. As the result, the orthogonality conditions transform to Mr = 0, 
where M is a positive and hence non-degenerate matrix. This implies r = 0. 
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